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SUMMARY 

An  airborne  survey  of  earth  electrical  resistivity,  computed  from  the  complex  till  ot  the 
electric  vector  ol  a VLI  (very  low  frequency)  surface  electromagnetic  wave,  has  been 
studied.  The  survey  was  analyzed  lor  the  elfects  of  topography  \s  resistivity,  altitude,  and 
wavetilt  phase  and  amplitude.  These  factors  can  seriously  affect  the  geologic  resolution  of 
this  type  of  survey.  Topographic  relief  al feels  at  least  one  electi  ic  field  component,  flight 
altitude  often  varies  over  relief,  and  phase,  which  cannot  he  separated  from  amplitude  b\ 
an  airborne  antenna  system,  is  also  quite  variable  when  the  earth’s  resistivity  is  stratified  or 
when  displacement  currents  are  dominant. 

A mountainous  area  in  northern  Maine  was  selected  for  study.  The  geology  and  topogra- 
phs of  the  area  had  been  previously  mapped  and  the  area  was  within  the  surface  wave  range 
ol  the  transmitter  used  (NAA,  Cutler,  Maine  17.8  kHz).  A standard  wavetilt  survey  was 
then  conducted  at  I 50-m  mean  flight  altitude.  Resistivits  effects  upon  the  wavetilt  field 
components  were  separated  from  those  of  topography  by  increasing  the  altitude  of  specific 
flight  paths.  I he  effect  of  altitude  upon  geologic  resolution  and  the  mean  measured 
resistivity  level  were  studied  by  repeating  part  ol  the  standard  survey  at  a higher  altitude. 

The  influence  of  phase  and  amplitude  upon  the  measured  airborne  wavetilt  was  analyzed 
using  the  results  of  a ground  survey  of  the  total,  complex  surface  impedance  ol  the  same 
radiowaves. 

The  results  provided  the  following  information.  The  igneous  lock  areas  generally  con- 
tained more  high  resistivity  anomalies  than  did  the  slate  areas.  Without  exception,  all  these 
anomalies  fell  on  mountain  flanks.  The  low  resistivity  areas  (ponds)  were  strongly  affected 
by  the  neighboring  high  resistivity  areas.  At  altitudes  above  300  m,  large  resistivity  responses 
remained  in  the  horizontal  electric  field  but  not  in  the  vertical  electric  field.  A topographic 
reevaluation,  based  on  an  adjustment  ol  all  values  of  the  vertical  electric  field  to  the  strength 
measured  over  areas  ot  no  relief,  improved  the  geologic  resolution  of  the  profiles  and  ol  the 
contouring.  The  repe.  ted  standard  survey  at  300  m retained  most  ot  the  geologic  resolution 
of  the  I 50-m  standard  survey.  I lie  ground  values  of  apparent  resistivity  and  phase  showed 
that  the  slate  areas  were  of  comparable  resistivity  to  the  igneous  areas  but  ol  much  lower 
phase.  They  also  showed  that  mountain  ridges  were  comparable  to  flanks  in  airborne 
apparent  resistivity  when  evaluated  by  the  topographically  corrected  method. 

Some  of  the  conclusions  reached  were  as  follows: 

1 . The  response  ol  the  vertical  electric  field  to  changes  in  resistivity  is  negligible  com- 
pared with  the  response  of  the  horizontal  electric  field  to  these  same  changes. 

2.  Topography  primarily  affects  the  vertical  electric  field,  causing  high  resistivity  zones 
to  appear  concentrated  on  flanks. 

3.  The  phase  of  the  wavetilt  is  an  important  factor  for  distinguishing  high  resistivity 
areas  in  regions  of  near-surface  crystalline  bedrock. 

4.  Resistivity  anomalies  less  than  twice  the  survey  altitude  in  dimension  were  difficult 
to  detect. 

5.  The  quality  of  a survey  may  improve  when  conducted  at  an  altitude  at  which  topogra- 
phs does  not  affect  the  vertical  electric  lie-id  and  will  improve  at  standard  altitudes  when  all 
values  of  the  vertical  electric  field  arc  adjusted  to  the  measured  strength  over  areas  ol  no 
relief. 
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Subscripts 

a apparent,  airborne 
cal  calibrated  value  of  a quantity 
i in-phase  value  of  a quantity 
m number  of  sky  wave  hops 
meas  measured  value  of  a quantity 

o free  space  value  of  a quantity  unless  otherwise  indicated 
q quadrature  phase  value  of  a quantity 
r receive 

s surface  vdue  of  quantity 
t transmit 

x,  y,  i components  of  a vector  quantity 
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INVESTIGATION  OF  AN  AIRBORNE  RESISTIVITY  SURVEY 
CONDUCTED  AT  VERY  LOW  FREQUENCY 


Steven  A.  Arcone 


INTRODUCTION 

Background 

Geophysical  methods  are  often  used  for  the  subsur- 
face exploration  of  new  sources  of  raw  materials. 

These  methods  utilize  the  physical  properties  of  rocks, 
sediments  or  certain  minerals  as  a basis  for  the  differ- 
entiation or  identification  of  material.  To  evaluate 
these  properties,  the  response  of  a particular  system  is 
usually  compared  to  that  of  a simulated  model  of 
idealized  geometry  and  material  properties.  Since  ideal 
conditions  are  rarely  encountered,  however,  further 
investigation  of  present  techniques  in  nonideal  environ- 
ments and  development  of  new  techniques  are  needed. 

Geophysical  exploration  often  uses  the  physical 
property  of  electric.il  resistivity  to  differentiate  materi- 
al types.  Materials  important  to  the  construction 
industry,  such  as  sands,  gravels  and  certain  cry  stalline 
rocks,  can  exhibit  resistivities  one  to  two  orders  of 
magnitude  greater  than  those  of  the  more  common 
materials  with  which  they  may  be  found.  Sands  and 
gravels  are  usually  sought  in  areas  of  formerly  existing 
deltas  and  floodplains,  where  the  usual  geological  lay  - 
ering and  a lack  of  topographic  relief  often  conform  to 
the  geometric  idealizations  of  mathematical  models 
used  tor  data  comparison.  Cry  stalline  rocks,  however, 
are  usually  associated  with  nonideal  mountainous  ter- 
rain. I herelorc,  before  large  amounts  of  time  and 
money  arc  invested  in  the  exploration  for  crystalline 
rocks  using  a geophysical  technique  based  on  an  ideal- 
ized flat  earth  theory,  careful  investigation  should  be 
made  of  the  response  of  that  technique  to  the  topo- 
graphic conditions  expected.  In  this  report,  the  per- 
formance of  an  airborne  resistivity  exploration  tech- 
nique over  a topographically  complex  region  is  experi- 
mentally analyzed. 

lire  geophysical  technique  discussed  uses  radiowaves 
propagated  Irorn  distant  transmitters  in  regular  use. 

The  field  vectors  of  these  waves  over  homogeneous 


earth  arc  shown  in  Figure  la.  In  this  figure,  / x is  the 
electric  field  component  tangential  lo  (he  earth’s  surface 
and  oriented  toward  the  transmitter,  / t is  the  vertical 
electric  field  component,  and  //v  is  the  magnetic  field 
component,  also  tangential  to  the  earth's  surface.  Be- 
neath the  surface  the  large  refractive  index  ol  the  earth 
allows  only  /.  and  Hs  to  exist,  as  the  wave  propagates 
vertically  downward  while  attenuating  exponentially. 

1 he  quantities  related  to  the  resistivity  of  the  earth  are 
the  surface  impedance  /s  defined  as 


and  the  wavelilt  If  defined  as 


Generally  , both  these  quantities  convey  the  same  in- 
formation and  either  one  or  the  other  is  measured.  Since 
antennas  can  be  used  to  measure  these  fields,  the  con- 
tact resistance  problems  of  conventional  galvanic  tech- 
niques that  use  earth  contact  electrodes  are  avoided. 

Both  2.  and  If  are  complex  quantities  whose  ampli- 
tudes and  phases  contain  resistivity  information.  The 
mathematical  theory  of  plane  waves  incident  upon  a 
“flat”  earth  (i.e.  with  topographic  relief  dimensions  that 
are  a small  fraction  of  the  wavelength  used)  containing 
horizontal,  homogeneous  layers  (Cagniard  1953)  is  the 
basis  for  converting  or  If  into  resistivity  and,  ulti- 
mately, geological  information.  An  "apparent  resistivity" 
p usually  expressed  in  ohm-meters  and  derived  from 
either  | Zs|  2 or  | If  j ’ is  most  commonly  determined.  At 
any  particular  frequency  p , represents  some  sort  of 
average  for  the  layering  to  a depth  which  increases  with 
decreasing  frequency . 


Remote 

Transmitter 
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a.  Field  vectors  above  homogeneous  earth  shown  retracting  vertically  into  the 
ground. 
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h.  Major  VLF  field  changes  known  to  take  place  when 
local  changes  in  earth  conditions  are  encountered. 

Figure  1.  I.  lee  tromagnetic  field  vectors  of  a radiowave  propagating  over  the 
earth 's  surface. 


Ideally,  then,  many  frequencies  could  be  used  to 
analyze  a resistivity  stratification.  Practically,  how- 
ever, only  transmitters  operating  in  the  VLF  band 
(very  low  frequency,  10-30  kHz  or  10-30  km  wave- 
length) are  powerful  enough  to  allow  their  radiation 
fo  be  used  worldwide.  When  radiation  at  higher  fre- 
quencies is  available,  the  information  provided  is  usu- 
ally more  difficult  to  interpret  because  of  the  undesira- 
ble effect  of  displacement  cutrents.  Transmitters  at 
lower  frequencies  do  not  exist;  therefore,  natural 
sources  of  radiation  in  the  LLF  (txtremely  low  fre- 
quency, 1-10 ’ Hz)  and  in  the  ULI  (ultra  low  frequency, 
I0_  ,-I  Hz)  bands  are  used.  The  primary  sources  within 
this  range  are  atmospheric  electric  discharges  and 
diurnal  variations  in  the  earth’s  magnetic  field  caused 
by  solar  emission,  ionospheric  emission,  etc.  The  use 
ol  this  radiation  is  known  as  the  magnctotelluric  (MT) 
method  (Keller  and  I rischknecht  1966,  Telford  et  al. 
1976).  In  this  report,  attention  will  be  confined  to 
only  VLF. 


Observed  VLF  responses  to  local  changes  in  an  earth 
environment  are  depicted  in  Figure  lb.  The  incident 
wave  may  be  a surface  wave*  or  a sky  wave  at  grazing 
incidence.  is  primarily  affected  over  gradual  hori- 
zontal changes  in  earth  resistivity  (Hoekstra  et  al.  1974) 
so  as  to  maintain  continuity  ol  current  flow  in  the  x 
direction.  //v  is  primarily  affected  over  large  lateral 
discontinuities  such  as  highly  conductive  ore-bearing 
dikes  (Paal  1965,  1968)  in  which  strong  eddy  currents 
may  be  induced.  £,  is  primarily  affected  over  large 
topographic  disturbances  (Harrison  el  al.  1971),  an 
electrostatic  effect  related  to  the  mountain  geometry 
and  therefore  occurring  independently  of  the  value  of 
earth  resistivity.  It  thcrelorc  violates  the  important 
assumption  of  a flat  earth  used  in  the  mathematical 
theory  employed  for  resistivity  data  reduction. 

All  the  above  responses  arc  generally  oversimplified 
because  changes  in  one  field  component  must  be 

* Also  referred  to  as  a ground  ssavc  in  the  literature. 
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accompanied  by  changes  in  all  other  field  components 
as  predicted  by  Maxwell's  equations.  However,  to 
date,  no  studies  have  been  made  to  demonstrate  that 
topographic  effects  can  be  separated  from  resistivity 
effects  for  any  VLF  field  component.  Consequently, 
resistivity  is  presently  computed  from  the  ratio  of  Ls 
to  i . or  to  //v  as  measured  in-situ. 

Presently,  Zs  is  commonly  measured  at  ground  level 
because  vegetation  interferes  with  E r Commercially 
available  instruments  are  capable  of  measuring  ampli- 
tude and  phase,  and  they  employ  as  the  phase 
reference.  At  airborne  levels,  above  vegetative  influ- 
ence, only  the  imaginary  or  quadrature  value  of  the 
complex  wavetilt*  IT'q  = Exq/En  provides  meaningful 
information  (E { is  used  as  the  phase  reference).  The 
reasons  for  this  will  be  discussed  later,  but  it  is  im- 
portant to  realize  now  that  amplitude  | M-' | and  phase 
angle  0 of  ft'  cannot  be  separated  by  an  airborne 
system,  since  W is  related  to  both  quantities  by  the 
relation: 

M'  = | W’j  sin  0. 

Therefore,  in  an  airborne  survey,  it  is  entirely  plausible 
that  a very  large  amplitude  of  ft’ caused  by  a very  re- 
sistive material  may  remain  entirely  undetected  because 
of  a low  value  of  0,  such  as  occurs  when  resistivity 
rapidly  increases  with  depth  or  when  dielectric  proper- 
ties become  important.  Also,  during  an  airborne  sur- 
vey, an  aircraft  often  changes  Might  altitude  by  several 
hundred  meters  when  navigating  over  steep  relief.  This 
in  itself  must  alter  the  geologic  resolution  of  a survey. 

Some  compensation  for  the  phase  effect  is  presently 
made  in  the  data  processing  by  assuming  that  0 is  usually 
near  45°  (Barringer  1972,  1973),  as  with  homogeneous 
ground  where  If,  = IT  and  the  true  resistivity  can  be 
determined.  Unfortunately,  homogeneity  is  not  com- 
mon over  the  depths  capable  of  being  penetrated  by 
VLf  radiation.  Although  Keller  et  al.  (1970)  and 
Hoekstra  et  al.  (1974)  have  assumed  that  45°  is  prob- 
ably correct  over  exposed,  deep-rooted  bedrock  forma- 
tions because  such  geology  itself  is  reasonably  homo- 
geneous, Sellmann  et  al.  (197$)  have  reported  Zs  phase 
values  as  low  as  I 3 over  exposed  deep  crystalline  forma- 
tions. Should  this  phase  hold  title  lor  IT,  then  these 
observations  would  contradict  the  4$L  phase  assumption 
tor  such  formations. 

Many  airborne  resistivity  surveys  have  been  made  in 
North  America.  Over  areas  of  minor  topographic  relief, 
they  have  usually  been  successful  in  delineating  regions 
of  permafrost  (Hoekstra  et  al.  197$,  Geological  Survey 


• a jnil  i refer  lo  qu.iilr.Hurc  .mil  in-phase  components,  re- 
spectively. 
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of  Canada  1973)  and  in  detecting  intertill  gravel  deposits 
(Culley  et  al.  197$).  All  verifications  have  been  made 
by  ground  reconnaissance,  and  in  the  work  of  Hoekstra 
et  al.  (1974,  1975)  ground  measurements  of  / were 
successfully  related  to  the  airborne  measurements  ol 
If  although  some  of  the  finer  detail  of  the  ground 
readings  did  not  appear  in  the  airborne  results. 

Of  the  few  surveys  known  to  have  been  conducted 
over  regions  of  more  severe  topographic  relief,  none  has 
been  adequately  verified  by  ground  reconnaissance. 

Keller  et  al.  (1970)  investigated  granitic  formations  in 
Colorado,  but  did  not  have  VLf  equipment  available 
tor  ground  resistivity  verification.  Both  Sellmann  et  al. 
(1974)  and  Hoekstra  et  al.  (1975)  conducted  surveys 
over  mountainous  regions  in  interior  Alaska  in  an  at- 
tempt to  delineate  permafrost,  but  made  no  ground 
observations  to  verify  their  conclusions.  An  examina- 
tion of  the  original  data  of  Hoekstra  et  al.  (1975)  by 
the  author  revealed  the  same  topographic  sensitivity  of 
E/  as  did  the  observations  by  Harrison  et  al.  (1971 ) 
mentioned  above. 

In  summary,  large  deposits  of  needed  raw  materials 
often  lie  in  mountainous  regions.  Since  these  deposits 
arc  usually  of  exceptionally  high  electrical  resistivity, 
a VLF  airborne  wavetilt  survey  would  seem  ideally 
suited  for  discovering  them.  However,  the  results  of 
previous  investigations  have  shown  that  three  basic 
assumptions  used  in  VLf  data  processing  a Hal  earth, 
45°-phase  angles  ovei  crystalline  bedrock  geology , and 
a steady  flight  altitude  can  be  seriously  violated. 
Consequently,  further  investigation  is  needed  into  the 
performance  ol  the  VLI  survey  system  in  mountainous 
terrain. 

Objectives 

The  experimental  objectives  were  to  determine  the 
effects  of:  1 ) topography  vs  resistivity,  2)  flight  alti- 
tude,! 3)  amplitude  of  the  wavetilt,  and  4)  phase  ol  the 
wavetilt  upon  the  geologic  resolution  of  a standard  air 
borne  wavetilt  survey  conducted  at  very  low  frequency. 
Based  on  previous  experiments  bv  I larrison  et  al.  (1971), 
it  was  expected  that  topography  would  severely  affect 
Ef j,  but  it  had  not  been  shown  either  experimentally  or 
theoretically  whether  topography  would  affect  / xq  or 
whether  changes  in  resistivity  would  aflect  I /r  Nor  had 
the  separate  influences  of  amplitude  and  phase  upon  the 
airborne  wavetilt  been  investigated  for  an  actual  survey . 
The  effect  of  flight  altitude  also  had  not  been  investi- 
gated. 


♦ In  this  report  the  term  “flight  altitude  " denotes  elevation 
above  terrain  and  not  above  sea  level. 


Experimental  procedures 

In  the  experiments  performed,  a mountainous  area 
in  northern  Maine  was  selected  for  study . The  area 
was  chosen  for  the  following  reasons: 

1 ) the  topography  was  mapped 

2)  the  bedrock  geology  was  mapped 

3)  the  overburden  was  usually  less  than  I m with 
common  outcrops 

4)  the  site  was  within  the  surface  wave  range  of  the 
transmitter  used,  and  therefore  the  radiation  was  free 
of  ionospheric  interference 

5)  apparent  resistivity  and  phase  contrasts  were 
known  to  exist  at  VLf  . 

The  area  was  survey  ed  at  VI J using  standard  airborne 
wavetilt  and  data  processing  techniques.  The  survey 
results  were  then  analyzed  as  follows  to  meet  the 
objectives  previously  mentioned. 

Topography 

The  mean  altitude  ot  specific  flight  paths  that 
traversed  the  steepest  relief  was  varied  to  eliminate 
topographic  effects  upon  the  lield  components.  The 
responses  at  the  higher  altitudes  provided  information 
as  to  which  field  components  were  responsive  to  re- 
sistivity changes  and,  hence,  which  were  responsive  to 
topographic  changes. 

Flight  altitude 

The  effect  of  flight  altitude  upon  geologic  resolution 
and  mean  apparent  resistivity  level  was  studied  by  nar- 
rowing the  flightline  spacing  of  part  ol  the  standard 
survey  and  then  repeating  it  at  twice  its  mean  flight 
altitude.  The  two  survey  s were  then  contoured  and 
compared. 

Amplitude  and  phase  ot  the  wavetilt 

These  factors  were  investigated  simultaneously  by 
comparing  apparent  resistivities  obtained  from  a ground 
survey  ot  the  total  complex  surlace  impedance  to  those 
obtained  from  the  airborne  wavetilt.  It  was  expected 
that  Hs  would  generally  be  in  phase  with  F f making 
this  procedure  valid.  As  will  be  shown,  L and  II.  were 
generally  in  phase  and  was  unaffected  by  topogra- 
phy. 

The  information  gained  from  the  high-altitude 
flights  justified  a topographic  correction  scheme  which 
was  devised  and  applied  to  the  original  survey.  Some 
ot  the  corrections  were  verified  by  the  results  ol  the 
ground  survey. 


t’RINCII’LLS  OT  LLLCT KOMAf.NET  1C 
RESISTIVITY  SURVEYING 

Resistivity  of  earth  materials 

The  physical  property  ol  electrical  resistivity  is  the 
basis  of  electromagnetic  lesistivitv  suiveving  tor  sub- 
surface exploration.  Resistivity  value*  id  most  earth 
materials  usually  tail  between  10  and  10,000  ohm-m. 

The  most  important  factors  that  determine  the  re- 
sistivity of  unfrozen  crustal  materialsare 

1.  Intrinsic  resistivity  ol  the  conslituient  minerals. 

2.  Textural  properties  such  as  grain  size,  porosity 
and  permeability . 

3.  Water  content. 

4.  Ionic  concentration  in  tree  water. 

5.  Ionic  concentration  in  adsorbed  water. 

factor  I is  usually  insignificant  for  most  materials. 

Common  minerals  such  as  micas  and  quartz  have 
resistivities  greater  than  10111  ohm-m  (Parkomenkho 
1967).  Ore  minerals  (e.g.  metallic  oxides  and  sulfides) 
are  much  more  conductive  but  arc  volumetricallv 
insignificant  in  most  materials.  Clay  minerals  (e.g. 
montmorillonite  and  kaolinite)  warrant  the  most  im- 
portant mineralogical  consideration  when  discussing 
resistivity  factors.  The  large  amount  of  ions  usually 
adsorbed  on  the  crystal  lattice  of  clays  (Grim  1953) 
produces  a surface  conduction  effect  that  can  appreciably 
lower  the  resistivity  ol  a sediment  or  rock  (Ward  1967). 

Factors  2 through  5 commonly  account  for  the 
resistivity  values  observed  in  earth  materials.  Conduction 
is  mainly  electroly  tic  and  takes  place  along  continuous 
lilms  ot  water  adsorbed  on  grain  boundaries  and  also 
through  water  that  may  fill  the  pore  spaces  available 
(Parkomenhko  1967).  Since  conductivity  (the  inverse 
ot  resistivity)  is  proportional  to  the  number  of  ions 
present  and  their  mobility,  the  ionic  concentration  of 
the  water  and  the  number  of  conducting  paths,  as  deter- 
mined by  the  material's  textural  properties,  are  of  ut- 
most importance. 

In  I igure  2 the  expected  ranges  of  resistivity  are 
compiled  tor  most  earth  materials,  based  on  past  sur- 
vey sand  laboratory  measurements.  The  unconsolidated 
materials,  clay  and  marl*  through  loose  sands,  generally 
show  an  increase  in  resistivity  with  particle  size.  This 
may  be  expected  from  the  corresponding  decrease  in 
available  adsorbing  surface  area.  Moraine  or  well-mixed 
glacial  deposits  exhibit  a wide  range  ot  resistivities  due 
to  large  variations  in  material  ty  pe.  The  consolidated 
materials,  chalk  through  crystalline  rocks  (e.g.  granite, 
gneiss  or  schist),  generally  show  an  increase  in  resistivity 
with  decreasing  porosity  and  permeability.  Clay  content 

* See  glossurv  for  dctinil ions  of  geologic  terms  used. 
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Figure  2.  Resistivity  ranges  of  most  earth  materials  (after  Culley 
etal.  1975). 


may  also  bo  significant,  as  with  sandstones  (Ward  1967), 
and  the  large  amounts  of  iron  and  magnesium  rich 
minerals  in  basaltic  igneous  rocks  can  make  them  less 
resistive  than  granitic  igneous  types.  In  all  cases,  the 
effects  of  weathering,  iointing  or  fracture  may  signifi- 
cantly lower  the  resistivity  of  the  fresh  rock  due  to 
increased  permeability,  formation  of  clay  minerals,  etc. 

Radiowave  propagation  within  the 
earth-ionosphere  waveguide 

Electromagnetic  resistivity  surveying  utilizes  radio- 
waves propagating  within  the  earth-ionosphere  wave- 
guide. These  waves  propagate  by  two  important  modes 
which  begin  to  form  at  a distance  of  about  one  wave- 
length from  the  antenna.  Both  modes  will  be  consider- 
ed because  their  relative  ranges  can  aftect  the  planning 
of  a resistivity  survey. 

The  two  propagation  modes,  as  radiated  by  a verti- 
cal electric  dipole  antenna  situated  at  the  earth’s  sur- 
lacc,  are  illustrated  in  Figure  3.  The  "surface  wave” 
is  so  named  because  its  field  strength  decays  with 
height  above  the  earth’s  surface.  The  "sky  wave”  is 
actually  the  sum  of  many  modes  formed  between  the 
ionosphere  and  the  earth  and  exists  at  all  frequencies 
below  about  108  Hz.  Both  modes  arc  tiansversc 
magnetic  and  their  field  vectors  at  the  earth’s  surface 
may  be  represented  as  in  the  figure.  £ refers  to  the 
electric  field  components,  //  to  the  magnetic  field 
component  and  x,  y,  and  / to  a local  right-hand  Cartesian 
coordinate  system.  / is  perpendicular  to  the  earth's  sur- 
face and  x and  y are  tangential,  a and  y , described  be- 
low, refer  to  the  incidence  angles  of  the  first  order  sky 
wave  mode  upon  the  ionosphere  and  earth,  respectively. 


As  has  been  experimentally  shown  in  a number  of  in- 
vestigations reported  by  Hollingworth  (1926)  and  by 
Alpert  (1963),  both  modes  do  not  necessarily  combine 
in  phase  at  the  earth's  surface.  However,  depending  on 
distance  from  the  transmitter,  usually  one  or  the  other 
mode  predominates. 

The  surface  wave  has  less  range  than  the  sky  wave 
and  is  relatively  more  important  at  distances  of  less 
than  400  km  at  VLl  to  the  transmitter.  At  a distance 
cl  such  that 

d 50/f1^  miles  (Iordan  and  Balmain  1968)  (1) 

or 

r 80 If' I*  km 

where  f is  in  megahertz,  Sommerfeld  (1909),  Weyl 
(1919)  and  Norton  (1936)  have  shown  that  surface 
wave  field  strength  is  determined  by  the  expression 

F,  - £,  Fjd  (2) 

where  L , is  the  vertical  electric  field  strength  in  V/m, 

is  the  initial  field  strength,  and  £ is  known  as  the 
Sommerfeld  attenuation  function.  £ depends  on 
flic  frequency,  distance,  and  electrical  properties  of  the 
earth,  and  graphs  are  available  for  finding  its  absolute 
value  for  both  electric  and  magnetic  dipole  type  antennas 
(Kraichman  1970).  At  VLl  frequencies,  the  higher 
power  of  existing  transmitters  can  allow  the  surface  wave 
to  extend  to  over  1000  km  over  sea  water  (Watt  1967). 

Beyond  about  800  km  from  the  transmitter  the  sks 
wave  predominates  over  resistive  earth  (Watt  1967). 
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Figure  3.  Two  radiowave  propagation  modes  in  the  earth-ionosphere 
waveguide  (see  text  for  explanation  of  symbols).  The  electromagnetic 
field  components  are  characteristic  of  radiation  from  a vertical  electric 
dipole. 


Call 

Location 

Frequency  (kH?) 

NAA 

Cutlet , Maine 

17.8 

NL.K 

|im  Creek,  Washington 

18.6 

NSS 

Annapolis,  Maryland 

21.1 

NBA 

Balboa,  Canal  Zone,  Panama 

24.0 

GBR 

Rugby , 1 ngland 

16.0 

UMS 

Corks , USSR 

17.15 

NW( 

Northwest  Cape,  Austialia 

23.3  and  IS. 5 

NI’M 

1 ualualei,  Hawaii 

2 M 

figured.  Coverage  for  well-known  Vff  stations.  Circles  enclose  approximate 
areas  w here  held  strengths  exceed  100  g \ ",- 
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The  sky  wave  may  then  be  represented  as  the  sum  ol 
incident  and  reflected  traveling  plane  waves  which  com- 
bine to  produce  many  standing  wave  modes  between 
the  ionosphere  and  the  earth.  Following  Watt  (1967), 
the  total  sky  wave  field  at  VLI  may  be  given  as  the  sum 
ol  all  possible  rays  which  may  undergo  one  or  more 
"hops"  as  illustrated  in  f igure  3.  f or  a particular  hop, 
the  vertical  electric  field  /: /m  as  transmitted  and  re- 
ceived by  vertical  dipoles  at  the  earth's  surface  may  be 
written  as 

L,m  - ^ FiKmCmF, 

\ cos  V ^1  + — 'j  1 e‘/(’,I-w/M  (3) 

where  / - initial  field  strength 

m = number  of  hops 

q = angle  ol  incidence  upon  the  ionosphere 
measured  Irom  the  vertical 
V - angle  ol  incidence  upon  the  earth  measured 
from  the  earth  tangential 
G,  = cos-  y/  = antenna  directivity  gain  factor 
/ , - transmitting  antenna  field  strength  pattern 
modification  produced  by  the  proximity 
of  a finitely  conducting  earth 
/ , - same  as  / , but  lor  the  receiving  antenna 
R - ionospheric  reflection  coefficient 
Cm  - ionospheric  locusing  lactor 
Ail  ■ difference  in  length  between  the  skv  mode 
path  and  gieat  circle  aic  between  the  trans- 
mitter and  receivei 
A = Irec  space  wavelength. 

At  distances  ol  less  than  2000  km  the  single  hop  sky 
wave  is  the  piedominant  mode.  The  ctleclive  range  of 
most  of  the  world's  VL.I  transmitters  is  shown  in 
f igure  4. 

Rjdiowavc  interaction  with  the  earth 

General  cons iiJerul ions 

Not  all  radiowave  energy  ol  both  the  sky  and  surface 
waves  is  confined  to  the  earth-ionosphere  waveguide.  A 
small  traction  is  retracted  below  the  earth’s  surface,  as 
represented  in  Figure  3.  In  the  earth  the  energy  is  partly 
absorbed  and  partly  scattered  back  to  the  surface  from 
the  subsurface  layers  of  rocks  and  sediments.  The  depth 
of  penetration  increases  with  decreasing  frequency,  allow- 
ing values  between  If)'  and  10*’  11/  to  be  ideally  suited 
lor  exploring  the  composition  of  the  first  lew  bundled 
meters  of  the  earth's  crust.  In  this  frequency  range, 


resistivity  is  the  predominant  electrical  property  that 
allows  differentiation  ol  material  type,  although  per- 
mittivity may  also  be  important  as  will  be  shown  be- 
low. 

These  electrical  properties  can  be  related  to  the 
wavetilt  h defined  as 

U hi 

and  to  the  surface  impedance  Zs  defined  as 


where  £x  and  L / arc  in  V/m,  II  ^ is  in  ampere-turns/meter 
and  Zs  in  ohms.  These  definitions  automatically  apply 
to  the  surface  wave,  but  the  sky  wave’s  grazing  angle  of 
incidence  y must  be  a;  proaching  0°  for  these  quantities 
to  be  useful.  E s need  not  be  in  phase  with  either  E/  or 
H. . , making  IF  and  Zs  complex.  Frischknecht  (1973) 
has  stated  that  the  amplitudes  and  phases  of  these  ratios 
for  the  sky  wave  at  VLF  are  not  affected  by  ionospheric 
events  or  by  the  total  path  of  propagation  from  the 
transmitter,  but  that  they  are  primarily  determined  by 
the  local  electrical  properties  at  the  point  of  measure- 
ment. Although  this  has  not  yet  been  verified  by  diurnal 
observations,  the  proximity  ol  y to  0 at  large  distances, 
the  extensive  range  of  existing  transmitters,  and  the  deep 
earth  penetration  make  VLF  best  suited  for  resistivity 
surveys. 

The  basic  definitions  used  in  electromagnetic  re- 
sistivity surveying  are  adopted  from  the  simply  layered 
Hat  earth  model  shown  in  Figure  5.  Often  this  model 
is  not  exactly  applicable,  but  it  has  proven  to  be  a useful 
guideline  for  data  interpretation  (Palacky  el  al.  1975, 

N ibelani  and  Rankin  1969).  Below,  some  ol  the  mathe- 
matical formulations  will  be  reviewed.  In  this  model  we 
find  ft  Zs/Z0  where  the  f ree  space  impedance  Z0  is 
defined  as 

Z(j  = \/q0/e0  - 377  ohms 

with  free  space  magnetic  permeability,  4;r  x 10' 
henry  /meter  and  the  free  space  electric  permittivity 
c0  = 8.85  - Hr12  farad/meter. 

I his  result  is  obtained  because  as  iy  approaches  0 

F,H<s  -*o-  <6> 

Therefore,  for  these  cases,  a measure  ol  either  IF  or  Zs 
will  provide  the  same  information.  The  derivation  ol 
the  following  formulas  may  be  lound  in  Wait  (1962). 


for  deriving  the  definitions  of  apparent  resis- 
tivity ( t -layer  case)  and  also  used  in  modern  which  equals  the  true  resistivity  when  a -1.  F or  such 

data  interpretation  techniques  (3-layer  case).  cases 


The  magnetic  permeability  p is  usually  taken 
as  equal  to  that  of  free  space  for  each  layer. 


(13) 


Homogeneous  earth  and  definitions 
of  apparent  resistivity 

When  the  electrical  properties  arc  uniform  with 
depth,  the  surlace  impedance  is  given  by 


where 


(XI 


and  the  phase  angle  0 of  either  Zs  or  It'  is  45°. 

For  resistivity  values  below  10’  ohm-m,  | H'j  ■ 0.12 
at  VLF.  Therefore,  serious  leakage  of  C/  into  the  hori- 
zontal antenna  that  measures  f\  will  occur  when  the 
antenna  is  mounted  on  an  unstable  platform  as  in  an  air- 
craft. In  this  case,  since  W and  Z are  complex  quanti 
ties,  only  the  component  ol  that  is  in  quadrature 
phase  with  either  Et  or  H can  be  used.  Using  the  wave- 
tilt  as  an  example, 

U = j H'j  (cosO  + /'  sinb),  (14) 

so  that  Ft'  the  quadrature  wavetilt  value  is 


is  the  refractive  index  of  the  earth.  In  this  expression, 


Ft’q  = j Ft’ j sin  0. 


(15) 


cu  = 2itf(f=  frequency  in  FI/) 

/ = V-T 

k = relative  permittivity 
p - resistivity  in  ohm-m. 


The  quantity  a,  defined  as 


p cukcq 


(9) 


is  a measure  of  the  relative  strengths  of  displacement 
and  conduction  currents  and  is  usually  much  greater 
than  unity  at  VLI  and  at  resistivities  less  than  10'* 
ohm-m. 

f or  all  practical  cases/;  I and  eq  7 may  be  simple 
tied  to 


(10) 


A measurement  of  only  Ft ' cannot  distinguish  the 
separate  effects  of  amplitude  and  phase  and  a phase 
angle  ot  45°  must  be  arbitrarily  assumed.  The  formula 
for  converting  It  q to  apparent  resistivity,  defined  as  pq 
in  this  case,  is 


(16) 


which  gives  the  correct  value  of  resistivity  lor  a homo- 
geneous earth.  When  phases  arc  different  than  45°,  pq 
is  related  to  pd  In  the  formula: 


Pq  = 2pj  sin20.  (17) 

As  was  shown  in  I igure  3,  the  refracted  waves 
propagate  nearly  vertically  into  the  ground.  The 
propagation  factor  lor  the  fields  within  the  earth 
is  exp  (/  \/-iu)p0lp  /),  which  can  be  rewritten  as 
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r iyure  6.  Skin  depth  of  radiowaves  as  a function  of  frequency  for  various  values 
of  resistivity. 
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Figure  7.  Apparent  resistivity  p , and  phase  Oat  JO  k / //  for  various  homogeneous 
earth  models  of  resistivity  p and  relative  permittivity  k.  both  pd  and  0 da  reuse  as 
k ini  reuses.  g>  is  plotted  in  equiphase  < ontours  and  approai  lies  as  as  k approac  lies 
I.  At  a value  of  k <0  and  p 30,000  ohm-m,  an  apparent  resistivity  of  18,000 
ohm-m  and  a phase  of  about  JO  will  he  determined  tram  the  / s or  II  method. 


Figure  ft'.  Quadrature  value  of  apparent  resistivity  pc|  (computed  from  eg  / 7)  at  20 
kHz  for  various  homogeneous  earth  models  of  resistivity  p and  relative  permittivity 
k.  At  a value  of  k aO,  an  actual  resistivity  of  30,000  ohm-m  will  be  measured  by 
an  airborne  system  us  an  apparent  resistivity  of  about  2,500  ohm-m. 


exp  [(/+ 1 ) x'capy/lp  z j and  is  valid  for  a I . A 
measure  of  the  depth  of  penetration  of  the  waves  is 
then  given  by  the  skin  depth  formula: 

ft  = sJlpiZapQ  , (18) 

which  is  the  depth  at  which  the  radiation  will  attenuate 
to  e'1  of  its  original  value  at  the  surface,  ft  is  plotted 
m I igure  6 as  a function  ol  frequency  tor  vat  ious  re- 
sistivities. Any  change  in  the  resistive  properties  ol 
the  earth  below  this  depth  will  not  produce  any  signifi- 
cant change  in  either  the  amplitude  or  phase  of  Zs  or 
It  When  changes  take  place  above  this  depth,  then 
both  amplitude  and  phase  will  be  affected. 

When  o I or  less,  n of  cq  8 tends  toward  a real 
quantity  suppressing  both  phase  and  amplitude  of 
and  ft  . this  is  illustrated  irr  I igure  7 lor  a frequency 
ol  20  ktl/  in  the  VI. I range.  I or  lixed  values  ol  re- 
sistivity, the  decreases  in  apparent  resistivity  and  phase- 
arc  plotted  as  a function  of  the  increasing  dielectric 
constant  k Ioi  various  homogeneous  earth  models  of 
resistivilv  p.  The  most  seriously  affected  quantity  is 
p of  eq  17  which  is  plotted  against  the  same  variables 
in  I igure  8.  I or  a particular  value  of  p,  pt]  decreases 
more  rapidly  with  increasing  k than  doespa  ol  I igure 
7 due  to  the  added  phase  factor.  This  is  significant  lor 
airborne  sy  stems  which  only  measure  p . Values  of 
k 20  arc  not  uncommon  for  many  crystalline  rock 
types  (Paikomenkho  1907). 


I 'ertically  stratified  homogeneous  layers 

When  earth  consists  of  two  or  more  layers  of  differ- 
ent thicknesses  and  electrical  properties.  If  and  Zs 
will  change  in  both  magnitude  and  phase,  depending 
upon  all  the  earth  parameters  involved.  Wait  (1962) 
presents  a simple  method  for  generating  the  formulas 
for  any  number  of  layers,  and  present  data  processing 
generally  uses  a three-layer  model  (Palacky  1975),  but 
only  a two-layer  case  will  be  considered  here.  This  case 
is  sufficient  to  demonstrate  some  ol  the  more  important 
effects. 

When  the  upper  layer  is  characterized  by  a resistivity 
p | and  depth  rt  and  the  semiinfinite  lower  layer  has 
resistivity  pi  and  a ■ 1 lor  both  layers,  then  the  sur- 
face impedance  is  found  from  the  formula: 


^s 


y/tu>p0p  i 


s/pslp  | + / tan  k , t , 

I + iy/pslp  | tan  A-  j r | 


(19) 


where  the  wave  propagation  numbei  k j is  delined  as 
k , = V-  'copo/p,  . 


An  apparent  resistivity  is  then  delined  in  the  manner  ol 
eq  ! I and  1 2 as 


(20) 
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a.  Apparent  resistivity. 


I iqure  9.  / he  apparent  resistivity  (a)  and  phase  (h)  of  two-layer 
earth  models  as  a function  of  layer  depth  at  a frequency  of  20 
kll/. 


which  is  correct  when  becomes  infinitely  small.  In 
l inure  9 are  plotted  the  variations  in  apparent  resistivity 
and  phase  lor  two-layer  earth  as  a function  ol  layer 
thickness  at  20  kHz.  When  the  upper  layer  is  relatively 
more  resistive,  the  phase  angle  lies  between  45  and  90  ; 
when  the  upper  layer  is  less  resistive,  the  phase  angle 
lies  between  0 and  45  . This  phase  behavior  can  be 
generalized  lor  any  resistive  structure  where  the 


resistivity  either  progressively  decreases  or  increases  with 
depth. 

When  tt  of  eq  9 is  comparable  to  unity  lor  the  upper 
layer,  and  p,  p , , the  phase  angles  are  strongly  pushed 
toward  90°  while  decreasing  the  apparent  resistivity. 
When  a is  comparable  to  unity  tor  the  lower  layer  and 
Pi  p | , the  phases  will  strongly  decrease  toward  0r 
while  also  decreasing  the  apparent  resistivity. 
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Irregular  surface  and  subsurface  conditions 

Extreme  mathematical  difficulties  arise  in  handling 
irregular  conditions  because  of  the  mixed  geometries 
that  are  always  involved.  Numerical  techniques  ((ones 
and  Price  1970,  Coggon  1971 , Swift  1971 ) have  been 
applied  for  studying  the  effect  of  lateral  inhomogeneities 
(e.g.  ore  bearing  dikes)  upon  plane  waves,  but  these  have 
been  primarily  at  magnetotelluric  frequencies  below  105 
Hz.  At  VLF  Lytle  (1975)  and  Hughes  and  Wait  (1975) 
have  investigated  cases  involving  a flat  earth  with  non- 
uniform  layering.  In  particular,  Hughes  and  Wait  have 
investigated  TM  waves  incident  at  values  of  \j/  approach- 
ing 0°.  They  demonstrate  that  both  the  phase  and 
amplitude  of  Zs  and  W qualitatively  correlate  with 
changes  in  subsurface  layer  dimensions  that  occur  over 
distances  far  shorter  than  a free  space  or  in-situ  wave- 
length. However,  their  results  also  show  that  the  phases 
of  these  two  quantities  are  not  necessarily  the  same, 
nor  arc  | lf|  and  |ZS|  always  related  by  Z0.  Wait 
(1971)  discusses  similar  “recovery”  effects  tor  ground 
waves  and  reviews  the  literature  on  this  subject. 

Layering  interpretations  of  simultaneous  measure- 
ments of  pl(  at  two  or  more  frequencies  by  an  airborne 
wavetilt  system*  have  been  compared  to  the  known 
geology  over  flat  regions.  Culley  et  al.  (1975)  success- 
fully correlated  the  interpretation  ol  pq  profiles  taken 
at  frequencies  in  the  BCB,+  LI  **  and  VLf  bands  with 
drill  hole  results  in  layers  of  till  and  gravel  in  Saskatche- 
wan, Canada.  Hoekstra  et  al.  (1975)  demonstrated 
the  sensitivity  of  both  VLF  and  LL  to  changes  in 
permafrost  conditions  in  several  study  areas  near 
1 airbanks,  Alaska.  An  examination  of  Hoekstra’s 
iriginal  data  by  this  author  revealed  that  at  VL1  7 1 
was  tar  more  stable  than  L\  over  apparent  resistivities 
ranging  between  about  100  and  4000  ohm-m. 

Ground-based  observations  have  been  limited  by  the 
availability  and  portability  of  suitable  equipment. 
Pioneering  work  in  wavetilt  measurements  of  ground 
electrical  properties  was  done  by  Lliassen  (1957),  King 
(1968)  and  Blomquist  (1970).  However,  extensive  use 
of  the  wavetilt  method  for  ground  level  geophysical 
surveying  is  limited  by  the  interference  of  vegetation 
with  £ Consequently,  the  surface  impedance  method, 
when  used  with  a system  developed  by  Collett  and 
Becker  (1967),  has  made  extensive  ground  surveying 
possible  at  VLf.  In  Hoekstra’s  work  mentioned  above, 
airborne  measurements  of  pq  derived  Irom  If  q were 
successfully  compared  to  limited  ground  measurements 
of  pq  derived  from  Zs  in  areas  of  marginal  relief  by  using 


* t or  the  earlier  development  ol  this  system  see  t risshknecht 
( 1971),  Barringer  (1072),  McNeill  et  al.  ( 197  1). 

) BC.B:  Broadcast  Band,  SOO-1600  kHz. 

**  LF:  Low  frequency,  J00-400  kHz. 


eq  17.  Since  the  airborne  measurements  showed  that 
Et  was  stable,  this  implied  that  /■/  was  also  stable. 

Thus,  the  relation  of  eq  6 was  probably  maintained,  al- 
though plane  wave  sky  mode  theory  (the  transmitter 
was  over  2400  km  away)  predicts  changes  in  all  com- 
ponents when  p changes.  Hockstra's  work  therefore 
demonstrated  that  either  ground  or  airborne  resistivity 
mapping  in  regions  of  little  relief  is  usually  determined 
primarily  by  fluctuations  in 

Observations  over  topographically  complex  areas 
have  been  limited  to  VLL.  Harrison  et  al.  (1971),  using 
antennas  mounted  below  a helicopter,  observed  that  the 
vertical  electrical  field  (Station  NAA,  Cutler,  Maine)  can 
intensify  over  mountain  ridges  by  more  than  6 db  (this 
effect  was  illustrated  in  Fig.  1 ),  but  that  the  effect  dis- 
appears at  flight  altitudes  comparable  to  the  mountainous 
relief  itself.  Harrison  theorized  that  this  was  an  electro- 
static effect  (the  mountain  size  is  a small  fraction  of  a 
VLF  wavelength)  and  reproduced  this  same  behavior  by 
superimposing  the  field  of  an  induced,  hemispherical 
dipole  model  of  the  mountains  observed  upon  the  field 
of  a VLF  ground  wave.  Harrison  used  the  same  electrical 
parameters  for  all  of  his  ground  models.  Sellmann  et  al. 
(1974)  observed  that  high  values  of  airborne  resistivity 
were  most  prevalent  in  valley  regions  of  a survey  near 
Fairbanks,  Alaska,  and  reasoned  that  this  resulted  from 
the  presence  of  permafrost.  A subsequent  examination 
of  these  data  revealed  that  the  vertical  field  (Station 
NLK,  jim  Creek,  Washington)  was  enhanced  over  ridge 
tops  as  with  Harrison’s  observations,  thereby  suppressing 
the  apparent  resistivity  of  the  mountainous  areas.  These 
data  also  revealed  that  the  horizontal  field  £v  was  more 
consistent  with  known  distributions  of  permafrost. 
Sellmann  et  al.  (1975),  using  the  surface  impedance 
method  over  a variety  of  near-surface  igneous  and 
mctamorphic  rocks  in  areas  of  variable  relief  in  northern 
Maine,  reported  a wide  range  in  both  apparent  resistivity 
and  phase.  At  apparent  resistivity  values  between  4000 
and  30,000  ohm-m,  the  phases  ranged  from  1 3°  to  58  . 
Phase  values  as  low  as  I 3°  can  seriously  affect  the  ability 
of  an  airborne  system  to  detect  such  high  resistivit 
values,  should  these  same  phases  hold  true  for  the  wave- 
tilt. 


INSTRUMENTATION 

Airborne 

The  airborne  E-Phaseft  wavetilt  system  was  used  for 
the  airborne  surveys.  In  this  system  illustrated  in  F igure 
10,  horizontal  and  vertical  electric  dipole  antennas  arc 


tt  t -Phase  (Barringer  Research  Ltd.,  Toronto,  Canada)  is  a 
patented  trade  name. 
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Figure  10.  Front  view  of  an  aircraft  system  of  crossed  dipoles  for  measuring  the  wave- 
tilt  components  f x and  / r The  aircraft  flies  orthogonally  to  the  transmitter  direction. 


mounted  on  a nose  cone  extending  from  a STOL 
(short  take-off  and  landing)  aircraft.  The  horizontal 
antenna  is  mounted  orthogonally  to  the  long  axis  of 
the  nose  cone.  The  aircraft  must  then  fly  perpendicu- 
larly to  the  transmitter  direction  for  maximum  coupl- 
ing with  E . This  allows  the  more  intense  Et  to  couple 
with  the  horizontal  antenna  when  the  aircralt  rolls  so 
that  only  the  total  E/  and  the  component  of  £x  in 
quadrature  phase  with  C { can  be  measured  success- 
fully to  determine  a wavetilt  value  H'q.  This  value  is 
then  determined  along  several  equispaced  tlightpaths 
over  the  survey  area. 


Calibration 

A wavetilt  calibration  is  performed  before  and  after 
each  flight  by  rotating  the  horizontal  antenna  90°  to 
the  transmitter  direction  for  minimum  coupling  with 
£fx.  This  antenna  is  then  tilted  20  to  the  vertical  to 
pick  up  a known  percentage  of  E r The  fcxq  channel 
is  delayed  90°  to  read  f xl  (in-phase  component)  and 
the  signal  thus  recorded  is  proportional  to  E/  sin  20\ 
The  wavetilt  IF  that  this  represents  is 


ft  = 


E,  sin  20° 

t, 


0.34. 


The  RF  gains  of  the  £xq  and  F channels  are  then 
set  to  give  sufficient  deflection  ol  a pen  trace  on  an 
analog  (chart)  recorder.  In-flight  changes  of  gain  in 
either  channel  are  noted  on  the  chart  recorder.  All 
data  are  also  recorded  in  digital  form  on  magnetic 
tape. 

After  a survey  H q is  computed  from  the  formula 

W QJi  < x ^vai 

q K'CJ,  sin  7 ^opv 


where  ft 


^xq-cab  ^ / -cal 

= E xq-mcas/fz-meas-  lhl'  Of  the  in-flight 

measured  fields 


7 = the  angle  between  the  average  true  azimuth 
of  a flight  line  and  the  true  azimuth  from 
the  center  of  the  survey  area  to  the  trans- 
mitter (for  maximum  coupling  7 = 90'  ) 

4ca]  = attenuator  setting  during  calibration 
A mCas  = attenuator  setting  during  the  survey. 

The  strengths  of  the  electric  field  quantities  are  re- 
corded in  volts  on  the  magnetic  tape  and  in  millimeters 
on  the  analog  chart  recorder. 

Flightline  spacing  and  attitude 

The  choice  of  a flightiine  spacing  depends  on  naviga- 
tional, scientific  and  economic  factors.  I lightlines 
closer  than  100  m to  each  other  arc  difficult  to  navigate, 
especially  over  topographic  rcliel  where  flight  altitude 
must  often  increase.  Usually  the  spacing  is  based  on 
the  field  strength  of  a vertical  dipole  which,  in  the  near 
field,  is  proportional  to  the  reciprocal  of  the  cube  of 
the  distance  from  the  antenna.  At  an  altitude  of  h m, 
the  intercept  of  the  -10  db  field  strength  contour 
(relative  to  the  ground  point  directly  below  the  aircraft) 
with  the  ground  is  a circle  of  approximately  2 h m 
diameter.  2h  is  then  used  for  the  flightline  spacing. 

Over  flat  areas/7  is  minimally  75  m.  For  particularly 
large  surveys,  the  spacing  must  often  be  widened  to  save 
on  expense.  In  this  case,  h is  then  increased  to  about 
half  of  the  flightline  spacing. 

Flightpath  recovery 

The  flightpaths  arc  recovered  with  on-board  pho- 
tography and  an  altimeter.  A fiducial  marking  system 
coupled  to  a clock  drive  exposes  the  fiducial  count 
onto  the  flightpath  camera  film  while  simultaneously 
recording  a fiducial  on  both  the  analog  chart  and 
magnetic  tape.  A 35-mm  flight  path  camera,  manu- 
factured by  I light  Research  and  equipped  with  a wide- 
angle  lens,  exposes  a continuous  film  strip.  The  flight 
number  and  time  of  dav  are  entered  manually  on  the 
analog  chart. 
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1 ho  aircraf  t height  above  the  terrain  surface  is 
determined  by  a Bon/er  TRN  70  radio  altimeter.  The 
output  of  this  instrument  is  recorded  on  both  the 
analog  and  digital  recorders. 

Data  processing 

The  first  step  in  the  data  processing  is  a smoothing 
ot  the  vertical  electric  field  with  a 1-Hz  low  pass  filter. 
Then,  assuming  an  arbitrary  wavetilt  phase  angle  of 
45  , If  is  converted  to  apparent  resistivity  pqj  using 
eq  16. 

pqa  is  also  tillered,  again  with  an  upper  cutoff  of 
I ffz.  Therefore,  at  an  average  flight  speed  ot  160 
km/h,  fluctuations  in  resistivity  occurring  approxi- 
mately every  45  m are  filtered  out. 

The  values  of  pqd  and  their  corresponding  fiducial 
marks  are  printed  out  in  digital  form.  The  pqj  values 
are  then  placed  upon  the  reconstructed  flight  lines 
that  are  superimposed  upon  an  aerial  photomosaic 
reconstructed  from  a composite  of  the  on-board 
photography.  A final  presentation  may  then  be  made 
in  the  form  of  contour  maps  (contouring  done  by 
hand)  or  as  resistivity  profiles  along  eacii  flightline. 

Ground 

Technique 

The  commercially  available  Geonics  TM16R*  was 
used  for  the  ground  surveys.  This  instrument  can 
measure  both  amplitude  and  phase  of  the  surface 
impedance.  In  Figure  1 1 the  instrument  is  depicted 
in  its  proper  orientation  for  measurement  with  respect 
to  the  appropriate  electromagnetic  field  quantities  of 
a ground  wave.  Coil  A is  first  used  to  castablish  the 
transmitter  direction.  Coil  A is  then  shut  off  and  the 
ratio  of  fv,  measured  between  two  probes  spread  10 
m apart,  and  //v,  measured  with  coil  B,  is  determined 
by  adjusting  the  amplitude  and  phase  dials  until  an 
inaudible  null  is  attained.  In  the  figure,  coil  A has 
been  tilted  from  the  horizontal  to  the  vertical  direction 
and  is  no  longer  coupled  to  //  . 

Tuning 

A variety  of  printed  circuits,  each  pretuned  to  a 
particular  V IT  frequency  in  present  use,  are  available 
with  the  system,  these  circuits  employ  phase-lock 
amplifiers  to  maintain  tuning. 

Null  detection 

The  achievement  of  a sharp  null  lor  an  accurate 
readout  depends  on  the  signal-to-noise  ratio  encounter- 
ed when  tuning  for  amplitude  and  phase  values.  For 

‘ Manufactured  by  t iconics  L id  . Toronto,  Canada. 
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Figure  / 7.  Schematic  of  the  Geonics  EMI6R  surface 
impedance  meter.  The  instrument  is  shown  in  its 
measurement  orientation  with  respect  to  the  radio- 
wave field  vectors  £ and  H . 

the  sky  wave  mode  this  ratio  depends  on  ionospheric 
conditions  as  well  as  the  total  path  of  propagation. 

I or  the  ground  wave  mode,  this  ratio  primarily  de- 
pends upon  distance  from  the  transmitter.  At  dis- 
tances within  the  ground  wave  range  (•  800  km  over 
earth)  a readout  accuracy  of  ±5  % for  amplitude  and 
±0.5%  for  phase  is  usually  obtained.  This  was  true  for 
the  surveys  performed  in  northern  Maine  which  were 
slightly  more  than  300  km  from  the  transmitter  used. 

Contact  resistance 

The  input  impedance  of  the  ground  probes  is  !08 
ohms,  0.5  x 10"'-  farad.  Therefore,  over  exposed  rock, 
a reliable  measurement  can  be  made  without  having 
the  contact  resistance  of  the  probes  influence  the 
measurement.  In  any  case,  the  probes  need  only  be 
laid  upon  the  ground  surface. 

Calibration 

The  EMI6R  is  factory  calibrated  by  applying  a 
known  voltage  between  the  electric  field  probes  and 
placing  coil  B within  a larger  coil  that  supplies  a 
magnetic  field  of  known  strength.  The  phase  dis- 
crepancy between  the  applied  voltage  and  the  magnetic 
field  is  adjusted  and  the  instrument  is  calibrated  for  each 
of  the  VLF  frequencies  desired.  The  calibration  con- 
siders that  the  hypothetical  electric  field  is  the  applied 
voltage  divided  by  the  measurement  probe  spacing  of 
10  m. 

The  instrument  used  was  bought  new  within  one 
month  o.  . e surveys  and  was  calibrated  prior  to  ship- 
ping. 
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Data  reduction 

The  amplitude  dial  of  the  EM16R  reads  directly  in 
ohm-meters  ol  apparent  resistivity.  II  needed,  the 
amplitude  of  the  surface  impedance  in  ohms  may  be 
recovered  by  inverting eq  II.  The  resistivity  scale  is 
logarithmic  and  ranges  from  0 to  30,001)  ohm-m.  The 
phase  scale  is  linear  and  ranges  between  O'  and  90  . 

This  phase  range  is  therefore  based  strictly  on  the 
plane-wave  model  described  in  Principles  of  Electro- 
magnetic Resistivity  Surveying  and  does  not  allow  for 
phase  discrepancies  between  t and  //s  outside  olHhis 
range. 

I or  the  comparison  ol  ground  to  airborne  resistivity 
measurements,  the  I M16R  data  can  be  converted  to 
an  apparent  resistivity  based  onl\  on  the  quadrature 
value  ol  the  surface  impedance  and  an  assumed  phase 
angle  c\  of  45  . I he  conversion  was  given  by  cq  I 7. 

METHODS 

Experimental  procedure' 

A topographically  and  geologically  well-mapped  site 
in  northern  Maine  was  selected  tor  studs.  Over  this 
site,  / 1 and  / ol  the  surlace  mode  radiation  from  a 
nearby  V 1. 1 transmitter  were  recorded  along  several 
equispaced  flightpaths  at  a mean  altitude  ol  I 50  m. 
These  field  components  were  then  compared  using 
standard  data  processing  techniques  to  produce  an 


apparent  resistivity  contour  map  of  the  site.  The  map 
was  then  compared  to  the  geology  of  the  region. 

I rom  this  standard  survey,  two  particular  flight 
paths  were  chosen  and  repeated  at  several  altitudes. 
Along  these  flights  berth  Et  and  /.  vq  were  measured 
and  compared  to  the  topography  and  geology  below. 
The  increase  in  altitude  attenuated  the  topographic 
effects  and  allowed  an  examination  of  only  the  resis- 
tivity information  contained  in  both  field  components. 
The  information  thus  gained  was  later  used  for  re- 
evaluating the  standard  survey  with  topographic  effects 
removed. 

Part  of  the  standard  survey  was  then  repeated  at 
twice  its  mean  altitude  to  study  the  effect  of  altitude 
upon  two-dimensional  geologic  resolution  and  mean 
apparent  resistivity  levels.  This  repeated  survey  used 
a lesser  flightline  spacing  than  the  standard  survey. 
Therefore,  a few  additional  flightlines  were  added  to 
the  standard  survey  so  that  both  surveys  had  an  equal 
number  of  flightlincs  over  the  same  area. 

An  extensive  ground  survey  of  determined  the 
individual  effects  of  phase  and  amplitude  upon  the 
standard  survey  results.  Traverses  were  made  over  all 
major  rock  types  present.  In  a few  cases,  the  ground 
readings  obtained  were  used  to  verify  the  topographic 
corrections  applied  to  the  airborne  readings. 

In  I able  I is  presented  a summary  ol  the  specifica- 
tions tor  all  the  airborne  and  ground  studies. 


Table  I.  Study  specifications. 


Mean  flight 

Mean  flight  path 

Mean  flight 

S tud\ 

altitude 

Vo.  of 

length 

spat  mg 

no. 

/ )t  st  ription  unJ  purpose 

(m) 

flights 

(km) 

(km) 

1 

1 ovs  altitude  standard  survey 
geologic  delineations  using 
standard  techniques. 

150 

13 

10 

0.4 

> 

High  altitude  series  A topo- 

90 

1 

10 

graphic  and  resistive  effects  upon 

150’ 

1 

10 

1 f and  /-  x of  study  1 . 

*00! 

1 

10 

600 

1 

10 

High  altitudes  series  B same 

90 

1 

10 

purpose  as  scries  A. 

150’ 

1 

10 

300  • 

1 

10 

600 

1 

10 

3 

Low  and  high  altitude  surveys 

1 50*  * 

10 

10 

0.16 

altitude  effects  on  geologic 
delineations  of  study  1. 

300 

10 

10 

0.16 

4 

t»round  survey  amplitude  and 
phase  effects  upon  geologic 
delineations  of  study  1 . 

0 

• Adapted  from  study  1. 

* Adapted  from  study  3. 

•*  Includes  6 flightlincs  of  study  1. 
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Figure  1 3.  Topographic  contours  and  elevations  of  the  survey  area.  Elevations  are  in  feet. 


Test  site  Deboullie  Mountain  area  and  is  currently  a state  park. 

The  area  is  accessible  only  by  numerous  logging  roads 
I oc  at  ion  which  are  passable  by  truck  between  about  1 )une  and 

The  area  chosen  for  the  experimental  survey  (see  1 November.  In  nearby  I rcnchville,  Maine,  there  is 

cover  photograph)  is  shown  on  the  map  in  Figure  I 2.  an  airport  which  was  used  for  the  mobilization  of  the 

It  is  located  approximately  19  km  southeast  of  the  con-  flight  crew, 
fluence  of  the  St.  )ohn  and  Allagash  rivers  in  the 


1962). 


Topography 

The  test  site  is  situated  in  a small  mountainous  belt 
that  continues  in  a north-south  direction  and  lies  at  a 
mean  elevation  between  1000  and  I 500  It  (300-450  m). 
The  topographic  contours,  simplified  from  the  U.S. 
Geological  Survey  topographic  map  for  the  fish  River 
Lake  quadrangle,  are  shown  in  I igure  I 3.  The  area 
contains  several  large  ponds  which  dominate  the  low 
elevations.  The  high  elevations  arc  dominated  by 
Gardner,  Whitman,  Deboullie,  and  Black  mountains, 
all  of  which  rise  between  700  and  800  It  (210-240  m) 
above  the  ponds. 


The  relief  is  often  severe  with  many  elevations  falling 
off  quite  rapidly,  from  Deboullie  Peak  to  the  shore  of 
Deboullie  Pond,  the  elevation  decreases  by  800  ft  (240 
m)  over  a horizontal  distance  of  approximately  1600 
ft  (480  m)  and  parts  of  this  flank  are  rock  slides.  Along 
other  flanks,  as  at  the  southwest  shore  of  Galilee  Pond, 
these  falloffs  are  actually  sheer  clifls. 

Geology 

The  bedrock  geology  of  the  test  site,  an  igneous 
stock  surrounded  by  slate  (designated  Seboomook* 
slate),  is  shown  in  Figure  14  (after  Boone  1962).  The 
intrusive  stock  is  divided  into  a northern  section  of  pre- 
dominantly syenite  and  a southern  section  of  granodi- 
oritc.  There  are  a few  smaller  inclusions  of  slate  in  the 
syenite.  The  southerly  slopes  of  Deboullie  and  Black 
mountains  are  syenite  while  the  northerly  slopes  are 
slate.  The  western  slope  of  Gardner  Mountain  is  slate 
while  the  eastern  half  is  syenite  and  granodioritc.  The 
western  slope  of  Whitman  Mountain  is  granodioritc 
while  the  eastern  slope  is  slate.  The  dip  of  the  slaty 
cleavage  is  commonly  greater  than  70“  and  the  strike 
is  northeast  to  southwest.  Geologic  profiles  at  depth 
have  not  been  made. 

The  surf  icial  geology  of  the  region  is  generally  a 
glacial  till  (McKim  and  Merry  1975)  of  less  than  1 .0  m 
thickness.  The  lower  elevations,  where  sediment  ac- 
cumulation would  be  expected  to  be  greatest,  are 
dominated  by  the  pond  chains. 

VLT  electromagnetic  characteristics 

Station  NAA  (17.8  kllz)  located  300  km  to  the 
southeast  near  Cutlet,  Maine,  was  monitored  lor  all  sur- 
veys. The  transmitter  is  a top-loaded  vertical  electric 
monopole  with  an  effective  radiation  ol  890  kW  of  con- 
tinuous power.  This  transmitter,  as  arc  most  other  VLf 
transmitters,  is  in  regular  use  by  the  United  States  Navy 
for  submarine  communication.  Information  is  trans- 
mitted by  amplitude  modulation  allowing  the  phase 
lock  circuitry  of  the  receiving  equipment  used  in  the 
experiments  to  stay  tuned. 

NAA's  distance  falls  within  the  limit  (307  km  at  17.8 
kHz)  defined  by  eq  1 and  thus  allows  an  evaluation  of 
the  relative  contributions  ol  both  the  surface  and  sky 
modes  to  the  total  vertical  electric  field  strength.  This 
calculation  is  important  because  it  demonstrates  that 
the  VLf  signals  received  were  not  subject  to  ionospheric 
conditions.  From  eq  2 and  3,  the  formulas  for  the  rela- 
tive contributions  of  both  surface  and  sky  modes  in 
decibels  above  I qV/m  as  derived  in  Watt  (1967)  are 


* Larlv  Devonian  Arc  Seboomook  is  a regional  designation. 
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L.  [db,  I juV/m|  = 109.5  + 10  log/3. 

- 20  logc/  + 20  log  |F|  (21) 

E/  ^ | db,  1 juV/m  | = 103.5  + 10  log  P, 

- 20  log  (d  + Ac/)  + 10  log  G,  (a,  i^) 

+ 20  log  F,  + 20  log  |ftm|  + Cm 

+ 20  log  F,  + 20  log  (cos  v»  (22) 

where  d anJ  Pr  are  expressed  in  kilometers  and  kilo- 
watts, respectively. 

These  formulas  are  for  ground-based  transmitting 
and  receiving  antennas. 

Considering  only  a one  hop  sky  mode  transmission, 
an  ionospheric  height  of  75  km,  an  ionospheric  re- 
sistivity of  10>  ohm-m,*  and  a maximum  average 
ground  resistivity  of  500  ohm-m,  the  terms  in  these 
equations  evaluate  as  follows: 

f\  = 890  kW 
F-  1.0 
At/  - 36  km 
a = 63° 

V - 27° 

G,  --  0.79 


20  log  | Rm  | -32  db 

Cm  0.4  db. 


Substituting  these  values  into  cq  21  and  22,  the 
relative  vertical  electric  field  strengths  at  a distance  of 
300  km  from  station  NAA  are 

E - 89.5  db,  1 pV/m. 

/ surface 

E,  = 59.5  db,  1 pV/m. 

'sky 

Since  E has  been  measured  at  well  above  66  db  above 
1 pV/m  (U.S.  Naval  Research  Laboratory,  signal  level 
contours,  station  NAA),  these  calculations  ensured 
that  only  the  surface  wave  mode  was  present  and  that 
interference  between  the  two  modes  and  ionospheric 
effects  would  not  be  of  any  significance. 

In  Table  II  is  presented  a summary  of  previous  VLF 
investigations  in  the  area  (after  Sellmann  cl  al.  1 97 5) 
taken  at  1 7.8  kHz. 

Since  these  phase  angles  arc  considerably  below  45°, 
it  was  expected  that  both  phase  and  amplitude  would 
play  a significant  part  in  the  results  of  an  airborne  sur- 
vey. 

Survey  traverses 

Study  no.  1 standard  survey 

The  flightlines  of  this  survey  arc  shown  superim- 
posed on  the  topography  in  Figure  1 5.  A mean  alti- 
tude of  1 50  m was  prescribed  lor  the  entire  survey;  at 
this  altitude  300  m is  about  the  optimum  flightline 
separation.  However,  the  altimetry  of  Hoekstra’s 
surveys  (lloekstra  et  al.  1974)  had  revealed  that 
pilots  tend  to  fly  too  high  when  encountering  steep 
relief.  Therefore,  the  spacing  was  prescribed  at  400  m 
which  then  required  13  flightlines  to  cover  the  area. 


Table  II.  Summary  of  previous  investigations  at  17.8  kHz  in  the  vicinity  of  Deboullie 
Mountain,  Maine  (after  Sellmann  et  al.  1975). 


\<>.  at 


Rot  k typi 

/ oi  at  ion 

measurement'. 

PA  ( 10*  ohm-m) 

0 ( decrees) 

Syenite 

Deboullie  Mtn. 

3 

9-30 

1 3-23 

(•ranodiontc 

Gardner  Mtn. 

10 

1 2-30 

15-40 

Whitman  Mtn. 

14 

3.8-27 

22-3S 

Slate 

1 oguc  Pond 

3 

10-1  1 

22-29 

* I his  is  based  on  the  average  values  ot  t .4  x 10^  H/  lor  ion 
collision  frequency  and  SX10®/m-  lor  electron  density  (Wan 
1967). 


19 


F igurc  1 5.  I hqhl paths  ot  the  standard  survey  conducted  at  a mean  altitude  ol  1 50  meters. 


As  shown  in  f igure  1 5,  the  flightlines  covered  all  the 
major  rock  types,  topographic  features  and  major  ponds, 
I he  slight  changes  in  flight  direction  were  caused  by 
attempts  to  navigate  cquispaced  llight-paths  (the  pilot 
fixed  on  predetermined  landmarks).  The  mean  azimuth 
bearing  of  the  flightlines  y with  respect  to  the  trans- 


mitter direction  was  81°.  As  the  Dcboullic  Mountain 
area  is  wilderness,  no  interference  from  any  power  or 
telephone  lines  or  any  other  man-made  structures  was 
possible.  I his  and  the  following  airborne  surveys  were 
conducted  between  22  September  and  5 October  1975 


Figure  16.  Flight  paths  of  high  altitude  series  A and  B. 


Study  no.  2 high  altitude  series 

Flightpaths  6 and  9 of  Figure  I 5 were  repeated  be- 
cause of  their  horizontal  position  over  Gardner, 
Deboullie  and  Black  mountains  and  because  they 
traversed  the  slate,  granodiorite,  syenite  and  the  ponds. 
Three  additional  altitudes  were  added  to  each  of  these 
flightpaths  and  arc  shown  superimposed  upon  the 
geology  in  Figure  16.  The  mean  altitudes  flown  were 
90,  1 50  (of  the  standard  survey),  300  and  600  m,  as 
indicated.  It  was  expected  beforehand  that  by  600  m 
all  topographic  influence  upon  the  VLF  fields  would 
disappear  in  accordance  with  the  observations  of 
Harrison  ct  al.  ( 1971 ) discussed  previously.  The  two 
separate  flight  groups  are  distinguished  as  series  A and 
series  B. 


Study  no.  3 high  and  low  altitude 
standard  surveys 

A more  intensive  survey  of  160-m  mean  flightline 
spacing  and  300-m  mean  altitude  over  the  central 
section  of  the  test  site  was  chosen  for  this  study.  The 
central  portion  chosen  was  suspected  to  be  well  within 
a high  resistivity  zone  so  that  the  lower  resistivity 
zones  that  may  have  existed  to  the  northwest  and 
southeast  of  the  survey  limits  would  not  affect  the 
300-m  survey  any  more  than  they  would  the  150-m 
survey.  Four  flightlines  were  then  added  between 
flightlines  5 and  10  of  the  150-m  standard  survey  to 
allow  comparison.  The  flightlincs  of  both  these  sur- 
veys are  shown  in  f igure  17  superimposed  over  the 
geology.  The  figure  reveals  how  difficult  it  was  to 
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Figure  1 7.  Flight  paths  of  the  (a)  low  and  (b)  high  altitude  standard  surveys. 


navigate  equispaced  flightlincs.  Therefore,  it  must  be 
emphasized  that  the  prescribed  line  spacings  are  mean 
values. 

Study  no.  4 — ground  surveys 

The  ground  traverses  made  in  relation  to  the  to- 
pography of  the  area  arc  shown  in  Figure  I 8 accompa- 
nied by  a legend  of  the  rock  types  encountered.  Meas- 
urements of  the  surface  impedance  phase  and  amplitude 
were  made  every  60  m.  At  a few  measurement  stations 
the  electric  field  probes  of  the  (iconics  LMI6R  were 
rotated  90c  to  find  it  the  tangential  electric  field 


deviated  from  the  transmitter  direction.  The  results 
were  always  negative.  The  word  "tangential"  is  used 
rather  than  horizontal  because  rarely  were  the  electric 
field  probes  exactly  on  a horizontal  plane.  Usually 
they  were  tilted  a few  degrees  to  the  vertical,  since 
many  measurements  were  made  along  Hanks  that 
dipped  down  from  the  transmitter  direction.  The  ap- 
parent resistivity  and  phases  over  the  ponds  are  inferred 
from  previous  measurements  over  New  Lngland  ponds 
where  ice  makes  them  accessible  in  the  winter  (Parrot 
and  Fleming  1970). 
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I igure  19.  Apparent  resistivity  pqj  contours  superimposed  upon  the  geology.  Shaded  areas  are  syenite  and 
granodiorite.  Contour  values  are  given  in  thousands  of  ohm-meters;  mean  flight  altitude  is  ISO  m. 


RESULTS 

Study  no.  1 standard  survey 

The  results  of  the  standard  survey  arc  presented  in 
f igures  19,  20  and  21.  I igures  19  and  20  show  the 
resistivity  patterns  superimposed  on  the  geologic  and 
topographic  bases.  These  patterns  were  hand  contoured 
from  the  digitized  resistivity  values  computed  along 


each  flightline.  The  contours  arc  given  in  lOOOohm-m 
intervals.  Also  shown  in  Figure  19  is  the  geologic  base 
by  itself  to  facilitate  the  comparison. 

The  most  important  features  of  Figure  19  are  the 
more  intense,  higher  resistivity  anomalies  that  are 
mostly  associated  with  the  granodiorite  and  syenite. 
Over  the  granodiorite  there  arc  anomalies  of  5000, 
6000,  8000  and  I 2,000  ohm-m.  Over  the  syenite  there 


Figure  20.  Anomalous  high  resistivity  /ones  superimposed  upon  the  topography.  Resis- 
tivity contour  values  are  given  in  thousands  of  ohm-meters;  elevation  contours  are  in  feet. 


arc  anomalies  of  9000,  10,000  and  1 1 ,000  ohm-m.  The 
slate,  which  occupies  the  largest  area  of  the  survey,  is 
generally  dominated  hy  the  2000  and  3000  ohm-m 
contours,  but  there  are  three  anomalies  of  5000  ohm- 
m (just  above  left  center),  8000  ohm-m  (just  below 


right  center  actual  value  is  8900)  and  1 2,000  ohm-m 
(just  north  of  Black  Pond).  The  pond  areas  contain 
no  high  anomalies,  as  should  be  expected,  but  they  are 
obviously  influenced  by  the  large  anomalies  that  border 
their  shores.  Lake  water  usually  has  resistivity  values 
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Figure  21.  Normalized  distributions  and  means  of  the  digitized  standard 
airborne  resistivity  values  pqa  for  the  four  major  materials  of  the  survey 
area.  Flight  altitude  is  1 50  m. 


of  less  than  200  ohm-m.  However,  the  larger  ponds 
very  gradually  deepen  to  more  than  20  m so  that  values 
near  1000  or  2000  ohm-m  are  not  unreasonable  near 
the  shores. 

The  location  of  the  high  resistivity  anomalies  in 
relation  to  the  topography  is  shown  in  Figure  20. 
Without  exception,  the  points  of  highest  resistivity 
fall  on  mountain  flanks  (note  that  the  aircraft  did 
not  avoid  the  peaks  see  flightpaths  shown 
in  Figure  I 5).  The  reason  is  that  the  lower  eleva- 
tions are  dominated  by  the  lower  resistivity  pond 
chains,  while  over  the  higher  elevations  (ridges)  the 
reference  field  E / was  enhanced,  as  will  be  shown  in 
the  next  section.  It  should  also  be  noticed  that  there 
arc  many  flanks  located  over  all  rock  types,  especially 
in  the  igneous  /one  south  of  Gardner  and  Deboullie 
ponds  and  just  west  of  Black  Pond,  which  are  not  un- 
usually high  resistivity  areas.  Above  about  6000  ohm- 
m there  is  a definite  differentiation  among  mountain 
flanks,  even  within  one  rock  type,  and  a definite  dis- 
crimination against  mountain  ridges.  Therefore,  the 
actual  coverage  of,  say,  the  igneous  geology  could  not 
he  correctly  inferred  from  Figure  20. 

f igure  21  gives  the  normalized  distributions  and 
means  of  the  digitized  resistivity  values  for  each  ol  the 


four  major  material  types  of  the  area.  The  rock  ranges 
generally  fall  within  the  broad  range  specified  for 
crystalline  rock  in  Figure  2.  As  expected  from  the 
contours,  the  slate  is  measured  as  the  least  resistive 
rock,  with  64%  of  all  values  falling  between  1000  and 
3000  ohm-m  and  a mean  value  of  2169  ohm-m.  The 
syenite  is  measured  as  most  resistive,  giving  the  highest 
average  and  the  heaviest  distribution  between  5000  and 
10,000  ohm-m.  The  fact  that  most  of  the  higher  air- 
borne resistivities  arc  associated  with  the  intrusive 
geology  seems  to  indicate  that  the  intrusives  are  far 
more  resistive  than  the  slate.  On  the  contrary,  it  will 
be  shown  that  the  slate  is  actually  of  comparable  re- 
sistivity and  that  the  hidden  factors  of  topography  and 
wavetilt  phase  often  determine  whether  a topographic 
or  geologic  area  is  mapped  as  an  airborne  anomaly  or 
not. 

In  the  following  sections,  the  separate  influences 
of  topography,  altitude  and  the  amplitude  and  phase 
of  the  wavetilt  upon  this  survey  aie  experimentally 
analyzed. 

Study  no.  2 - topographic  and  resistive 
effects  upon  Txp  and  F t j 

In  this  study,  topographic  effects  upon  ffvq  and  L fi 
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/ igure  22.  L'ncalibrated  altimetry  (broken)  and  normalized  f f.  (solid)  analog  traces  along  the  fliyhtpaths  of  alti- 
tude series  A.  The  traces  are  not  exactly  coincidental  with  each  other  or  with  the  topographic  profile  because  ot 
changes  in  flight  speed.  Note  that  in  the  90- m trace  t /j  tends  to  return  to  nearly  the  same  level  during  periods 
ot  stable  altimetry  irrespective  of  the  geology  and  that  at  600  m the  disturbances  have  nearly  disappeared. 


are  separated  from  resistivity  effects  when  the  altitude 
of  two  particular  tlightpaths  is  increased  to  values 
comparable  (300  m)  and  above  (600  nt)  the  mountain- 
ous relief  itself.  It  is  assumed  that  any  topographic 
effects  upon  / are  attenuated  with  altitude  as  rapidly 
as  those  upon  t. n,  allowing  only  the  resistivity  responses 
of  the  two  field  components  to  remain  at  the  higher 
altitudes.  Ihe  traces  ol  the  individual  lield  components 
are  not  calibrated  in  terms  of  absolute  field  strength 
because  transmitter  power  and  gain  levels  varied  over 
the  several  days  that  the  flights  were  carried  out. 

I igures  22-25  show  the  normalized  analog  traces 
of  t/1  and  l.  vq  and  the  uncalibrated  altimetry  for  the 
flightpaths  in  I igure  16.  Ihe  topographic  and  geologic 
cross  sections  are  shown  at  the  bottom  of  the  figures. 
The  altimeter  traces  are  uncalibrated  but  terrain  clear- 
ance increases  upwards.  The  altitude  levels  given  are 


approximate  mean  values.  The  different  profiles  within 
a particular  series  are  not  all  exactly  coincidental  with 
each  other  or  with  the  topographic  profile  because  of 
changes  in  flight  speed  that  resulted  from  whether  a 
flight  was  from  southwest  to  northeast  or  vice-versa. 
Therefore,  for  each  trace  the  altimetry  must  be  used 
to  recognize  where  a particular  disturbance  takes  place. 
In  all  figures  the  topographic  vertical  scale  has  been 
magnified  by  a factor  slightly  greater  than  10. 

Figures  22  (series  A)  and  23  (series  B)  show'  the 
traces  ol  Eei.  The  normalization  values  /. 0 correspond 
to  the  more  stable  altimetry  over  Gardner  and  Togue 
ponds  (f  ig.  22  and  23).  At  90-m  and  1 50-m  mean 
altitudes  t7j  increases  above  all  the  higher  land  eleva- 
tions (this  is  seen  more  clearly  by  the  coincidence 
of  the  altimetry  changes  with  the  t./l  changes),  in 
these  cases  very  near  the  mountain  peaks.  Since  pqi 
is  inversely  proportional  to  the  square  of  f n.  these 
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Figure  23.  Uncalibrated  altimetry  (broken)  and  normalised  I /j  (solid)  analog  traces  along  the  flight- 
paths  of  altitude  series  li.  The  traces  are  not  exactly  coincidental  with  each  other  or  with  the  topo- 
graphic profile  because  of  changes  in  flight  speed.  Note  that  in  the  90-ni  trace  E f - tends  to  return  to 
nearly  the  same  level  during  periods  of  stable  altimetry,  irrespective  of  the  geology. 


traces  now  reveal  why  the  anomalies  ot  Figure  20  were 
confined  to  the  mountain  flanks.  At  the  300-m  level 
the  fluctuations  have  become  much  less  intense  and 
they  have  practically  disappeared  at  the  600-m  level  of 
series  A.  It  should  also  be  noticed  that  along  the  90-m 
traces,  f/j  always  tends  to  return  to  nearly  the  same 
level  over  the  flatter  portions  of  the  topographic  pro- 
file, irrespective  of  the  geology. 

These  observations  of  L {[  are  similar  to  those  of 
Harrison  et  al.  (1971 ) as  mentioned  above,  but  for  this 
Hudy  they  def  ine  the  600-m  altitude  level  as  the 
approximate  limit  of  topographic  or  resistive  influence 
upon  En  and,  therefore,  the  assumed  limit  of  topo- 
graphic influence  (if  any)  upon  £xq.  In  the  £xq  traces 
of  f igures  24  (scries  A)  and  25  (scries  B),  the  normaliza- 
tion values  were  again  taken  above  Gardner  (Fig.  24) 
and  Togue  (Fig.  25)  ponds.  As  with  the  t.  r[  traces,  all 
the  perturbations  have  attenuated  with  height.  Flow- 


cver,  at  600  m large  resistivity  perturbations  stili  re- 
main. They  are  much  larger  than  the  perturbations 

the  600-m  L n traces,  especially  those  of  I igure  22 
where  there  are  only  negligible  changes  above  the  same 
topographic  features. 

Since  all  the  above  field  strength  traces  arc  un- 
calibrated in  terms  of  amplifier  gain,  the  effect  of 
resistivity  upon  only  £ at  an  altitude  above  topo- 
graphic influence  is  best  illustrated  in  the  apparent 
resistivity  profiles  (solid  curves)  of  Figures  26  (series  A) 
and  27  (series  B)  because  they  account  for  the  calibra- 
tion of  both  field  levels.  At  the  600-m  altitude  ot  both 
flight  series  the  resistivity  values  may  be  seen  to  vary 
between  1000  and  7000  ohm-m  lor  series  A and  be- 
tween 1700  and  5000  ohm-m  lor  series  B.  Since  pqa 
is  proportional  to  (£xq/£/j)  ’,  £xq  varies  by  as  much 
as  a factor  of  2.78  as  compared  to  1 .05  fot  l n at  600 
m for  series  A.  For  series  B l xq  varies  by  as  much  as  a 
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Figure  24.  Uncalibrated  altimetry  (broken)  and  normalized  £xq  (solid)  analog  traces  along  the  flightpaths  of  altitude 
series  A.  Note  that  at  600-m  mean  flight  altitude  the  fluctuations  in  F . ..  are  much  more  intense  than  they  were  in 
Figure  22  for  F.  r over  the  same  topographic  features. 


factor  of  2.06  compared  to  1.20  for  FZI.  It  therefore 
seems  ev  ident  that  when  topographic  effects  are  elimi- 
nated only  Esq  is  essentially  responsive  to  resistivity 
changes. 

It  may  be  noticed  in  these  figures  that  the  lower  alti- 
tude (90-  and  I 50-m)  profiles  are  not  entirely  consistent 
with  the  high  altitude  (300  and  600  m)  profiles.  In 
Figure  26  the  Gardner  Mountain  anomaly  is  too  low  at 
90  and  1 50  m,  and  at  1 50  m the  Black  Mountain 
anomaly  is  relatively  too  high  compared  with  the  same 
anomaly  in  the  other  three  profiles.  In  Figure  27  the 
mean  values  at  90  and  1 50  m are  inconsistently  low, 
since  they  arc  so  much  closer  to  the  igneous  area.  How- 
ever, a significant  improvement  results  when  the  re- 
sistivities are  recomputed  (dashed  curves)  using  only 
the  reference  area  levels  F0  of  Figures  22  and  23  for 
I n . Now  the  syenite  anomalies  at  90  and  I 50  m over 
Gardner  Mountain  in  Tigure  26  are  more  consistent 
with  those  at  the  higher  altitudes,  with  each  other,  and 


also  with  the  other  syenite  anomalies  over  Deboullie 
and  Black  mountains.  The  Deboullie  Mountain  anomaly 
at  1 50  m is  now  the  largest  at  this  altitude  as  it  is  at  all 
other  altitudes.  In  Figure  27  the  mean  values  for  the 
different  altitudes  arc  all  now  more  comparable  and  the 
anomalies  to  the  right  of  Galilee  and  Deboullie  ponds 
have  been  enhanced  at  the  90-m  level  to  their  corrected 
values  at  the  150-m  level.  For  both  series  the  300-m 
and  600-m  profiles  have  changed  marginally  as  is  to  be 
expected. 

Study  no.  3 - general  effects  of  altitude 
upon  geologic  resolution 

The  altimeter  traces  of  Figures  22-25  revealed  that  it 
was  very  difficult  to  maintain  constant  terrain  clearance, 
although  the  flight  altitude  relative  to  sea  level  was  prob- 
ably much  more  constant.  The  kind  of  relief  encountered 
can  present  a serious  danger  to  an  aircraft  trying  to 
maintain  a low  survey  altitude  and  the  variability  in 
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Figure  25.  Uncalibrated  altimetry  (broken)  and  normalized  l (solid)  analog  traces  along  the  llight- 
paths  of  altitude  series  B.  Vote  that  at  600-m  mean  flight  altitude  the  fluctuations  in  are  much 
more  intense  than  they  were  in  Figure  23  for  E over  the  same  topographic  features. 


altimetry  may  well  affect  the  results.  The  profiles  of 
f igures  26  and  27  presented  sulficient  evidence  that 
surveys  over  this  kind  of  terrain  may  be  conducted  at 
higher  and  safer  altitudes  without  serious  loss  in 
resolution.  The  300-m  flights  of  those  figures  were 
adapted  from  this  studs  which  is  a more  generali/ed 
study  of  the  effects  of  altitude  upon  a standard  re- 
sistivity survey. 

The  flightlines  of  the  surveys  discussed  were  given 
in  I igure  17  I he  flightlinc  spacing  was  narrowed  to 

0.16  km  (mean  value)  to  miprove  resolution.  In 
f igure  2f>  are  compared  the  resistivity  contours  of 
the  two  surveys  f rom  these  contours  the  following 
observations  may  be  made 

I.  All  the  major  anomalies  of  the  lower  altitude 
repeat  at  the  higher  altitude  except  lot  the  one  between 
Gardner  and  Toguc  ponds.  This  9000  ohm-m  anomaly 


at  1 50  m reduces  to  tvy.o  smaller  5000  ohm-m  anomalies 
at  300  m. 

2.  Not  all  the  high  resistivity  anomalies  at  I 50  m 
decrease  in  intensity  at  300  m.  The  8000  ohm-m 
anomaly  between  Gardner  and  Deboullic  ponds  in- 
creases to  10,000  ohm-m  and  two  new  5000  ohm-m 
anomalies  appear  above  (northeast  of)  Black  Pond. 

3.  The  contours  over  the  ponds  are  generally  of 
higher  value  at  the  300-m  level. 

figure  29  gives  the  normalized  distributions  of  the 
digitized  resistivity  computations  lor  all  flightlines  of 
both  surveys.  The  normalized  distribution  ol  the 
300-m  survey  is  more  concentrated  about  the  mean 
value,  which  is  very  similar  tor  both  surveys,  but  has 
a much  greater  percentage  19%  vs  1 0%  ol  values 
between  4000  and  5000  ohm-m.  This  increase  reflects 
the  numerous  5000  ohm-m  anomalies  of  the  300-m 
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/ iqure  26.  Standard  (solid)  and  topographically  corrected  (broken)  pl|a  profiler  tor  alti- 
tude series  I.  p is  the  mean  ol  the  standard  and  p’  is  the  mean  ot  the  corrected  pro- 
files. 


contours.  The  30()-m  distribution  also  has  a much 
smaller  percentage  3%  vs  14%  of  values  below 
1000  ohm-m  which  reflects  a decrease  in  sensitivity 
to  the  low  resistivity  areas. 

The  results  of  both  surveys  are  generally  very 
similar,  demonstrating  that  most  resistivity  informa- 
tion measured  with  a standard  airborne  wavetilt  sys- 
tem over  topographic  relief  this  severe  is  preserved  to 
at  least  a 300-m  mean  flight  altitude.  Also,  for  the 
1 5 0-m  standard  survey,  the  many  local  changes  ex- 
perienced in  flight  altitude  did  not  significantly  affect 
the  resistivity  contours  to  be  derived  from  the  standard 
data  processing  method. 


Study  no.  4 - ground  level  study  of  phase 
and  amplitude  effects 

I igurcs  30  and  31  are  the  profiles  ol  pls  and  along 
each  traverse  of  Figure  1 8.  These  values  were  used  to 
profile  pl)b  (derived  from  eq  1 7)  in  Figures  32  and  33 
which  will  he  compared  to  the  airborne  results  in 
F igurc  19  or  20.  It  is  assumed  that  1. 1 and  II s arc 
generally  in  phase,  thereby  making  the  comparisons 
valid.  The  p()s  profile  values  will  be  explained  with 
reference  to  the  data  of  Figures  30  and  31. 

The  ground  traverses  are  described  below.  The 
parenthesized  numbers  after  the  traverse  letters  refer 
to  the  number  of  stations  where  measurements  were 
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Figure  27.  Standard  (solid)  and  topographically  corrected  (broken)  pqj  profiles  for  altitude  series  />'. 
p is  the  mean  of  the  standard  and  p ’ is  the  mean  of  the  corrected  profiles. 


made.  The  subscript  s refers  to  the  ground  (surf ace) 
value  of  a quantity 

1.  Traverse  A- A' (1 2).  This  is  predominantly  under 
a 3000  ohm-m  region  of  the  airborne  survey  just  above 
Togue  Pond.  pqs  = 1 560  ohm-m  not  including  the 
measurements  at  the  pond  shore.  The  three  measure- 
ments taken  at  the  pond  shore  were  very  low  at  350, 

240  and  1 80  ohm-m. 

2.  Traverse  B-B’  (1 1 ).  Ibis  is  under  the  5000  and 
6000  ohm-m  airborne  contours  north  of  Denny  Pond 
in  the  granodiorite  section.  p ls  = 4473  ohm-m,  <5S  = 
43.5°  and  pqs  = 3978  ohm-m. 

3.  Traverse  C-C’  (II).  This  is  under  the  1 2,000 
ohm-m  anomaly  west  of  Denny  Pond  in  the  granodiorite 
section.  One  value  of  pqs  exceeds  I 2,000  ohm-m  but 

p is  much  less  at  5360  ohm-m.  p LS  1 1 ,545  ohm-m 
which  is  very  close  to  the  peak  airborne  value,  but 

= 30.6°  which  explains  the  lower  value  for  pqs. 

4.  Traverse  D-D'  (63).  This  runs  the  length  of 
Togue  Pond.  It  is  generally  under  the  2000  and  3000 


ohm-m  airborne  contours  and  p = 1 825  ohm-m. 

5.  Traverse  t-l  ’ (41 ).  This  covers  all  three  rock 
types  as  it  passes  up  Gardner  Mountain  from  south  to 
north.  In  the  granodiorite  section  it  passes  under  an 
8000  ohm-m  airborne  anomaly . Over  the  granodiorite, 
pd5  = 9009  ohm-m  but  pqs  = 4374  ohm-m  because 

<5S  = 28.9°.  Over  the  syenite  section,  which  is  directly 
on  the  ridge  of  Gardner  Mountain,  the  traverse  is 
under  3000,  2000  and  1000  ohm-m  airborne  contours 
but  pqs  is  much  greater  at  6100  ohm-m. 

6.  Traverse  T-l  ' (46).  The  first  part  of  this  traverse 
is  over  the  slate  and  has  a very  low  pqs  = 1646  ohm-m 
because  0S  = 16.9°.  Over  the  granodiorite,  <>s  increases 
to  35.9°,  pJS  = 7086  ohm-m  and  pqs  = 4597  ohm-m. 
Over  this  traverse  the  airborne  contours  are  3000, 

4000  and  5000  ohm-m. 

7.  Traverse  G-G’  (13).  This  is  over  slate  and  goes  up 
Whitman  Mountain  directly  under  an  airborne  8000 
ohm-m  (actual  peak  value  is  8900)  anomaly.  pq>  = 

8503  ohm-m  and  <5S  is  a relatively  high  41°. 
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Figure  33.  pi)s  profiles  at  I 7.8  kHz  along  ground  traverses  F-f  through  /■  / 


8.  Traverse  HU’  (25).  This  is  all  over  syenite  and 
passes  directly  under  the  10,000  ohm-m  anomaly  be- 
tween Deboullic  and  Gardner  ponds.  As  shown  in 

I igure  31,  four  of  these  readings  are  indicated  at 
30,000  ohm-m,  hut  they  actually  exceeded  the  cali- 
brated range  of  the  instrument  and  could  not  be  deter- 
mined. <5S  was  exceptionally  low  along  this  traverse 
at  21.4  with  a standard  deviation  of  only  22%  for  all 
25  readings. 

9.  I raverse  I- 1’  (II).  I Ins  is  directly  across  the  twin 
peaks  of  Deboullie  Mountain  and  at  the  edge  of  the 
syenite.  The  airborne  contours  indicate  values  be- 
tween 1000  and  3000  ohm-m  but  p()s  is  much  higher 

at  4685  ohm-m.  i\  is  high  at  42". 


10.  Traverse  ) -)’ (40).  This  traverse  passes  from 
slate  to  syenite  and  stays  under  a 3000  ohm-m  contour 
1 his  corresponds  well  with  pqs  2480  ohm-m  tor  the 
slate  and  p |s  = 3492  ohm-m  for  the  syenite. 

Phase  angles  near  45  do  occur,  such  as  on  the  east 
flank  of  Whitman  Mountain  and  over  the  ridges  of 
Gardner  and  Deboullie  mountains.  Most  usually , how- 
ever, they  are  between  20  and  35  for  all  rock  types. 
At  only  I 5 of  the  270  ground  stations  is  phase  ever 
greater  than  45  , and  the  average  of  these  1 5 readings 
is  48.7°. 

In  some  cases  the  correlation  of  pqj  with  p ^ is 
excellent  as  with  traverses  G-G’  and  )-)  ’.  The  syenite 
section  of  traverse  L-L'  and  all  of  l-l’  are  directly  over 
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Ft igure  34.  Normalized  distributions  and  means  of  the  ground  values  pd 
and  0 


ridges  where  L y is  suppressed.  Correspondingly,  pl|s 
is  well  above  pc|d.  In  the  cases  of  traverses  C-C’  and 
H-H’  the  ground  readings  arc  much  lower  than  the 
highest  airborne  reading  made  over  the  same  area. 

Figure  34  gives  the  normalized  distributions  and 
means  of  the  ground  readings  organized  according  to 
rock  type.  Generally,  all  three  types  have  similar 
distributions,  with  67%  of  the  syenite  values,  75%  ot 
the  granodiorite  values  and  70%)  of  the  slate  values 
falling  between  I and  10,000  ohm-m.  The  syenite  has 
the  highest  followed  by  the  granodiorite.  The  slate 
is  almost  as  resistive  as  the  granodiorite  and  lar  more 
resistive  than  the  airborne  readings  indicate.  Since  only 
1 3 of  the  1 27  slate  readings  were  taken  under  high 
airborne  anomalies  (all  on  Whitman  Mountain),  there 
is  no  bias  in  the  slate  data  for  these  areas. 

The  phase  distributions  of  I igure  34  show  more 
discrimination  among  the  rock  types  than  do  the  pJV 
distributions.  The  syenite  has  two  definite  peaks, 
one  between  20  and  25  and  another  between  40° 


and  45°.  The  former  is  mostly  due  to  flank  readings 
and  the  latter  to  ridge  readings.  Only  19%  of  the 
syenite  readings  are  within  ± 5°  of  the  mean.  The 
granodiorite  has  a single  strong  peak  between  35  and 
40°,  and  55%  of  the  readings  are  within  ± 5°  ol  the 
mean.  The  slate  has  no  distinct  peaks  above  10%  and 
has  the  lowest  mean.  Forty-seven  percent  of  all  slate 
phases  fall  below  25  as  opposed  to  only  I 1.5%  for  the 
granodiorite,  which  has  a comparable  pAS  distribution. 

Figure  35  compares  the  distributions  of  pqs  and 
pqj.  For  all  rock  types  the  phase  has  significantly 
lowered  the  ground  resistivities  so  that  the  three  ground 
distributions  arc  all  comparable  to  those  of  the  airborne. 
No  values  now  appear  above  16,000  ohm-m  for  the 
ground  values.  In  all  cases  pqa  is  much  closer  to  pqs 
than  it  is  to  pas.  Eighty-one  percent  of  all  airborne 
slate  values  are  below  3000  ohm-m;  this  compares 
favorably  with  71%  for  the  ground  values.  Seventy- 
three  percent  of  all  airborne  granodiorite  values  are 
between  2000  and  6000  ohm-m,  again  comparing 
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Figure  35.  Comparison  between  the  normalized  distributions  and  means  of  the  standard  airborne 
and  ground  surveys  of  pt). 
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a.  Airborne  and  ground  values. 


Figure  36.  Comparison  of  airborne  contoured  values  and  ground  readings  of  tor  tra- 
verse F-F".  Dots  and  circles  are  actual  measurements. 


b.  Location  of  traverse. 
Figure  36  (cont 'd). 


favorably  with  7 1 .5%  of  the  ground  values.  For  the 
syenite,  61  of  all  values  fall  between  2000  and  6000 
ohnt-m  for  the  airborne  survey,  compared  to  55%  for 
the  ground  survey.  The  mean  ground  values  given 
tor  each  rock  type  ate  slightly  higher  than  the  airborne 
values  and  this  is  probably  due,  in  part,  to  the  topo- 
graphic enhancement  of  F t discussed  previously. 
Later,  when  topographic  effects  are  eliminated  and 
the  standard  survey  is  reevaluated,  some  ol  these  com- 
parisons between  the  means  will  improve. 

Over  regions  where  topography  did  not  seriously 
disturb  the  airborne  results,  specific  comparisons 
between  ground  and  airborne  values  may  be  made  to 
approximately  determine  the  necessary  si/e  ol  a resis- 
tivity anomaly  for  airborne  detection  by  contouring. 


In  Figures  36  and  37  are  compared  pqs  along  traverses 
F-F'  and  J-J  ’ with  the  overhead  1 50-m  contoured 
values  of  pqj.  Points  and  circles  on  the  curves  indicate 
actual  measurements.  In  Figure  36  an  intense  ground 
anomaly  (peak  value  - 14,600  ohm-m)  exists  for  a 
length  of  about  300  m and  is  interpreted  by  the  airborne 
contouring  as  a 5500-ohm-m  anomaly.  To  the  right, 
a lesser  ground  anomaly  (peak  value  = 8400  ohm-m) 
also  occurs  for  a length  ot  about  300  m but  is  barely 
detected  by  the  contouring.  In  Figure  37  one  large- 
ground  anomaly  (peak  value  9000  ohm-m)  occurs 
at  the  right  for  slightly  more  than  300  m and  is  dis- 
cerned by  the  contouring.  In  both  figures,  ground 
anomalies  of  shorter  length  are  hardly  seen  at  alti- 
tude. Apparently  300  m is  near  the  signif  icant  physical 


Distance  (km) 

a.  Airborne  and  ground  values. 


I igure  37.  Comparison  between  airborne  contoured  values  and  ground  readings 
of  Pq  for  traverse  j-f  . Dots  and  circles  are  actual  measurements. 


b.  Local  ion  of  traverse. 

figure  37  (cont d).  Comparison  between  airborne  contoured  values  and  ground 
readings  ot  pq  for  traverse  Dots  and  circles  are  actual  measurements. 


dimension  which  ground  anomalies  of  these  intensities 
(above  background  levels  of  3000  to  4000  uhrn-m) 
must  have  to  register  as  an  airborne  anomaly  a!  I 50- 
m altitude  and  approximately  400-m  flightline  spacing. 
This  approximate  si/e  corresponds  to  the  - 10  dh 
sensitivity  width  of  an  airborne  dipole  at  l 50-m  alti- 
tude. 


THE  STANDARD  SURVEY  REEVALUATED 

\ he  studies  of  flight  series  A and  B in  Results  re- 
vealed that  only  E n is  disturbed  by  topographic  relief 
while  £ contains  the  actual  resistivity  information. 

I he  behavior  of  E/t  supports  the  earlier  observations 
of  Harrison  (1971),  discussed  previously.  The  ground- 
air  comparisons  ot  pq  over  mountain  ridges  (traverses 
E-E’  and  1-0  demonstrated  that  pqd  is  suppressed  in 
these  settings  by  this  effect. 

Along  each  flightline  (I  ig.  I 5)  ol  the  standard  sur- 
vey, L n responds  to  the  topographic  relief  in  the  same- 
manner  of  I igures  22  and  23.  Over  each  mountain 
ridge,  L n rises.  Over  the  llattcr  sections,  such  as  the 
ponds,  all  the  E n levels  are  within  a few  percent  of 
each  other  as  was  evident  in  I igures  22  and  23.  there- 
fore, at  each  fiducial  mark,  a new  apparent  resistivity 
pqd  was  computed  by  adjusting  / /|  along  an  entire 
flightline  to  the  field  level  above  these  Hat  sections. 

The  formula  for  pqd  is 


f,qa  = Pl,a(f/i/£'o)2  (23) 

where  (23) 

Pqd  = (2/cjc0)  (E^/E,y 

as  before  and  E0  is  the  value  ol  E ■ over  the  flat  areas, 
such  as  the  ponds.  This  undisturbed  level  is  not  neces- 
sarily the  minimum  level.  Often  £d|  dips  slightly  below 
the  undisturbed  level  on  both  sides  of  the  peak  response 
to  a topographic  disturbance. 

I igurc  38  presents  this  reevaluation  ot  the  standard 
survey  . The  locations  of  the  high  resistiv  ity  patterns 
are  similar  to  those  of  I igure  19  but  several  features 
have  been  changed: 

1.  One  anomaly  ot  f igurc  19  is  considerably  more 
intense,  having  risen  from  10,000  to  17,000  ohm-m. 
This  lies  on  a flank  between  Deboullie  and  Gardner 
ponds  and  also  occurs  where  several  ground  readings 
could  not  be  determined  because  they  exceeded  30,000 
ohm-m. 

2.  The  anomaly  over  Gardner  Mountain  has  now 
become  three  anomalies  that  include  the  entire  moun- 
tain. Above  the  syenite  section,  the  pl|d  contours  have 
risen  from  1000-3000  ohm-m  in  I igure  19  to  6000  and 
7000  ohm-m  in  f' igurc  38.  This  agrees  better  w ith  the 
value  pqs  ot  5342  ohm-m  tor  eight  ground  readings 
along  this  ridge,  and  the  corrected  values  along  a section 
ol  flightline  no.  9 over  Gardner  Mountain  are  shown  in 
Figure  39. 
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/ igurc  38.  Apparent  resistivity  contours  of  the  topographically  corrected  standard  survey  superimposed  up- 
on the  geology.  Shaded  areas  are  syenite  and  granodiorite.  Contour  values  are  to  he  multiplied  h\  1,000. 
Mean  flight  altitude  = 150  m. 


3.  Between  the  twin  peaks  of  Deboullie  Mountain 
(northern  boundary  of  the  syenite  section),  the  air- 
borne contours  have  risen  from  1000-2000  ohm-m  in 
f igurc  19  to  3000-3000  ohm-m.  This  agrees  well  with 
the  value  4685  ohm-m  for  I I ground  readings 
along  this  ridge. 


4.  The  anomaly  over  Whitman  Mountain  in  the  slate 
section  has  marginally  decreased  from  8900  ohm-m  to 
8100  ohm-m.  This  is  also  in  good  agreement  with  the 
value  b(|s  r 8500  ohm-m  for  I 3 readings  across  the 
southern  flank. 


a.  Comparison  of  standard  and  corrected  pq  values  for  flightline  9. 


Figure  39.  Apparent  resistivities  pqa  and  pqa  computed  along  flight  line  no.  9 over  Gardner  Mountain 
and  topographic  and  geologic  profile  beneath  the  flightline.  The  eight  ground  samples  in  (b)  were 
made  linearly  in  a direction  transverse  to  the  flightline.  The  corrected  values  over  the  syenite  are  more 
consistent  with  the  ground  readings  and  with  other  airborne  anomalies  over  the  syenite  regions  of 
Deboullie  and  Black  mountains. 


a.  a 000  ohm-m  contour  corrected  survey.  b.  4000  ohm-m  standard. 

/ igure  40.  Comparison  between  the  igneous  areas  encompassed  by  the  a 000  ohm-m 
contours  of  the  topographically  corrected  survey  (a)  and  the  standard  survey  (b).  Un- 
shaded regions  cover  the  igneous  geology. 
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Figure  4 1.  Comparison  betw  een  the  normalized  distributions  and  means  of  the  topographically  corrected 
airborne  survey  of  p'  and  the  ground  survey  of  p 


5.  In  the  top  portion  ol  the  figure,  4000  and  5000 
ohm-m  contours  now  appear  as  they  did  in  the  300-m 
survey  of  f igure  28  but  not  in  the  1 50-m  survey  of 
Figures  1 9 or  28. 

6.  A new  anomaly  has  occurred  in  the  slate  section 
on  the  western  approach  to  Deboullie  Mountain.  No 
ground  readings  were  made  to  verify  this. 

A more  general  comparison  between  the  standard 
and  the  corrected  survey  is  made  in  I igure  40.  In  this 
figure,  the  areas  encompassed  by  the  4000  ohm-m 
contour  are  compared.  This  contour  value  was  chosen 
because  it  is  just  below pqs  for  both  igneous  rock  types 
but  well  above  pqs  ol  the  slate  values.  Therefore,  this 
contour  might  be  expected  to  best  define  the  limits  of 
the  igneous  geology.  As  shown  in  the  comparison,  the 
topographically  corrected  contours  contain  much  more 
ol  the  igneous  geology  in  all  sections  of  the  map.  Only 
at  the  very  top  of  the  figure  is  this  contour  overextend- 
ed in  relation  to  the  uncorrected  standard  survey. 

fhe  normalized  distributions  ol  resistivity  values  tor 
the  reevaluated  survey  are  given  in  f igure  41 . The  com 
parisons  are  similar  to  those  ot  f igure  35  where  similar 
comparisons  were  made  with  the  standard  airborne  re- 
sults. In  this  case,  however , the  mean  values  for  the 
slate  and  granodiorite  have  drawn  closer  to  those  ol 
the  ground  measurements,  although  pql  lor  syenite  has 
increased  to  over  6000  ohm-m. 


CONCLUSIONS  AND  FINAL  REMARKS 

Important  conclusions  to  be  drawn  from  these  VLF 
resistivity  studies  are  as  follows: 

1.  Above  ridges  and  peaks  Etl  was  always  enhanced 
in  the  manner  predicted  by  the  theories  of  Harrison 

et  al.  (1971)  lor  electrostatic  dipole  models  of  moun- 
tains. As  verified  by  ground  measurements,  this  effect 
reduced  the  apparent  resistivity  (as  measured  with 
standard  airborne  procedures)  above  these  locations 
and  caused  the  higher  resistivity  zones  to  be  concentrat- 
ed strictly  on  mountain  flanks. 

2.  The  high  altitude  flights  showed  that  resistivity 
predominantly  affects  IS  while  E,,  is  relatively  un- 
affected by  resistivity. 

3.  Most  resistivity  information  was  preserved  to  at 
least  3(X)-m  altitude  when  the  resistivity  was  computed 
from  the  actual  E and  Cxq  using  standard  procedures. 
Ibis  resulted  from  the  compensation  for  the  amplitude 
loss  in  £xq  by  the  amplitude  loss  in  £/i  over  mountain- 
ous regions. 

4.  Geologic  resolution  improved  when  the  airborne 
apparent  resistivity  was  recomputed  using  the  topo- 
graphically corrected  L n. 

5.  Wavetilt  and  surface  impedance  phase  were  gen- 
erally similar,  as  deduced  from  the  correlation  ol  ground 
and  airborne  pq  values.  This  implies  that  E ; and  Hs 
were  also  generally  in  phase. 

6.  The  phase  of  the  wavetilt  was  an  important  factor 
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for  distinguishing  high  resistivity  areas  in  regions  of 
near  surface,  crystalline  bedrock. 

I . Resistivity  anomalies  less  than  approximately 
twice  the  survey  altitude  in  dimension  were  difficult 
to  detect. 

Recordings  of  £x|  were  also  made  but  were  highly 
erratic.  This  behavior  was  expected  from  the  coupling 
of  L 1 1 into  the  horizontal  dipole  during  aircraft  roll 
but  was  also  probably  due  to  topographic  scattering. 
Topographic  scattering  is  the  primary  reason  why 
antenna  stabilization  schemes  designed  to  measure 
the  full  wavetilt  cannot  work:  the  in-phase  compo- 
nent of  the  wavetilt  is  more  than  likely  destroyed  be- 
fore it  can  be  measured  by  an  airborne  antenna.  The 
phase  of  surface  impedance  correlates  with  wavetiit 
phase  only  as  determined  by  £X[)  and  £XI  before 
scattering.  It  seems  probable,  however,  that  L x|  at 
ground  level  is  a meaningful  quantity  because  so  many 
phase  measurements  were  between  20  and  35°, 
which  is  plausible  lor  the  geology  encountered. 

The  observed  phase  angles  of  less  than  45°  are 
probably  due,  in  a general  sense,  to  a combination  of 
the  two  factors:  I ) a stratification  in  resistivity  such 
as  would  be  caused  by  a gradation  in  weathering  with 
depth,  and  2)  high  permittivity  values,  especially  for 
the  slate  for  which  values  as  high  as  34  have  been  re- 
ported by  Parkomcnkho  (1967).  The  latter  factor 
would  imply  that,  even  at  VLF,  crystalline  rock 
distinctions  could  have  actually  been  made  on  the 
basis  of  relative  permittivity.  It  permittivity  were  a 
factor,  then  some  of  these  results  might  he  different 
it  another  transmitter,  located  in  a direction  orthogonal 
to  the  direction  of  Cutler,  Maine,  had  been  monitored. 

I his  is  because  slate  can  be  anisotropic  due  to  its 
cleavage.  A study  using  orthogonal  transmitters  in 
such  an  area  might  then  define  surface  impedance  ten- 
sors from  which  additional  geologic  information  might 
possibly  be  uncovered. 

I he  theoretical  basis  ol  this  exploration  method 
demands  infinite  layers  ol  homogeneous  material  in 
order  to  correlate  the  lield  ratios  with  the  resistivity. 

At  the  area  studied,  the  largest  horizontal  dimension 
ol  the  igneous  region  was  less  than  half  of  the  wave- 
length (16.8  km  at  17.8  kHz)  monitored.  Therefore, 
it  might  seem  improbable  that  lield  measurements 
processed  on  the  basis  of  a theory  that  can  be  derived 
horn  geometric  optics  should  provide  results  consistent 
with  values  determined  by  more  conventional  field  and 
laboratory  techniques.  Nevertheless,  these  studies 
present  observational  evidence  that  ground  waves  can 
rapidly  recover*  (in  terms  of  fractions  of  a wavelength) 

*Ki  t,*fcncr  to  the  theofv  *'t  the  ''recoverv '*  i-ttrct  was  given 
rv  hr,  ij.ild'  surfjt e and  subsnrtdt » i -mdit  >>ns 


the  desired  relative  electromagnetic  field  strengths  and 
pha>es  when  passing  from  one  resistivity  condition  to 
another. 
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APPENDIX:  GLOSSARY  OF  GEOLOGIC  TERMS  USED 

(adapted  from  Glossary  of  Geology  by  the 
American  Geological  Institute,  Washington,  D.C.) 


Basalt 

Chalk 

Clay 

Cleavage 

Crystalline 

Dip 

/ erromagnesian 
r dilation 

Glassy 

Gneiss 

Granite 


A dark  to  medium-dark,  commonly 
volcanic,  igneous  rock. 

A soft,  pure,  fine-textured  usually 
white  to  light  gray  or  huff  limestone 
of  marine  origin,  consisting  mainly 
of  calcite. 

A loose,  earthy,  natural  sediment  or 
soft  rock  composed  primarily  of 
clay-si/e  (<  1/256  mm)  particles  and 
characterized  by  high  plasticity  and 
mostly  hydrous  aluminum  silicates. 

The  property  or  tendency  of  a rock 
to  split  along  secondary,  aligned 
fractures  or  other  closely  spaced, 
planar  structures  or  textures,  pro- 
duced by  deformation  or  metamor- 
phism. 

A descriptive  term  referring  to  a rock 
consisting  wholly  of  crystals  or  frag- 
ments of  crystals,  e.g.  an  igneous 
rock  lacking  glassy  material,  or  a 
metamorphic  rock. 

The  angle  that  a structural  surface, 
e.g.  a bedding  or  lault  plane,  makes 
with  the  horizontal,  measured  per- 
pendicular to  strike  of  the  structure. 

Containing  iron  and  magnesium. 

A general  term  for  a planar  arrange- 
ment of  textural  or  structural  features 
in  any  type  of  rock,  e.g.  cleavage  in 
slate  or  schistosity  in  schist. 

A descriptive  term  referring  usually 
to  the  texture  or  luster  of  usually  a 
volcanic  rock.  Glassy  properties  re- 
sult from  rapid  cooling  without 
distinct  crystallization. 

A foliated  metamorphic  rock  in 
which  bands  of  granular  minerals 
alternate  with  bands  of  I la ky  miner- 
als. 

A term  loosely  applied  to  any  light- 
colored,  coarse  grained  igneous  rock. 
Specifically,  it  refers  to  a definite 


composition  of  quartz,  feldspar  and  fer- 
romagnesian  minerals. 

Granodiorite  A lyp*-  of  granitic  rock  ol  specified  com- 
position of  quartz,  feldspar  and  lerro- 
magnesian  minerals. 

Igneous  A descriptive  term  referring  to  a rock  or 
mineral  that  solidilied  from  a molten 
(magmatic)  or  partly  molten  material. 

Intrusive  A descriptive  term  usually  referring  to 
an  igneous  rock  that  has  solidified  from 
a magma  emplaced  in  other  rock 
structures. 

Limestone  A sedimentary  rock  consisting  chiefly  of 
calcium  carbonate. 

Loam  A rich,  permeable  soil  of  similar  pro- 

portions of  clay,  silt,  sand  and  organic 
matter.  Somct  mes  called  topsoil. 

Marl  A term  loosely  applied  to  unconsolidated 

materials  consisting  chiefly  ol  clay  and 
calcium  carbonate. 

Metamorphic  A descriptive  term  applied  to  a rock 

derived  from  pre-existing  rock  by  min- 
cralogical,  chemical,  and  structural  changes 
in  response  primarily  to  marked  changes  in 
temperature  and  pressure. 

Permafrost  Any  material  in  which  a temperature  be- 
low 0°C  has  existed  continuously  for 
more  than  about  two  years.  The  defini- 
tion is  based  exclusively  on  temperature. 

Sandstone  A sedimentary  rock  of  sand  size  particles 
set  in  a silty  or  clayey  matrix  and  cement- 
ed usually  by  silica,  iron  oxide  or  calcium 
carbonate. 

Schist  A strongly  foliated  metamorphic  rock 

which  can  be  readily  split  into  thin 
flakes  or  slabs;  i.c.  it  demonstrates  the 
property  of  schistosity  which  is  a type 
of  cleavage. 

Slate  A fine-grained  metamorphic  rock  that 

cleaves  along  planes  independent  of  the 
original  bedding. 


Stuck  An  igneous  intrusion  less  than  about  40 
miles-  in  surface  exposure. 

Strike  I he  direction  or  trend  that  a structural  sur- 
face, e.g.  a bedding  or  fault  plane,  lakes  as  it 
intersects  the  horizontal. 

Syenite  An  igneous  rock  similar  to  granite  but  of  less 
quartz  content. 

lill  Unsorted,  unstratified  and  unconsolidated 

mixture  of  clay,  sand  and  gravel.  I ill  is  usually 
of  glacial  origin. 
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f or  conversion  of  SI  metric  units  to  U.S. /British 
customary  units  of  measurement  consult  A STM 
Standard  L380,  Metric  Practice  Guide,  published 
by  the  American  Society  for  Testing  and  Materials, 
1916  Race  St.,  Philadelphia,  Pa.  19/03. 
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